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NONHOLONOMIC LAGRANGIAN SYSTEMS ON LIE 
ALGEBROIDS 


JORGE CORTES, MANUEL DE LEON, JUAN C. MARRERO, AND EDUARDO MARTINEZ 


Abstract. This paper presents a geometric description on Lie algebroids of 
Lagrangian systems subject to nonholonomic constraints. The Lie algebroid 
framework provides a natural generalization of classical tangent bundle ge¬ 
ometry. We define the notion of nonholonomically constrained system, and 
characterize regularity conditions that guarantee that the dynamics of the 
system can be obtained as a suitable projection of the unconstrained dynam¬ 
ics. The proposed novel formalism provides new insights into the geometry 
of nonholonomic systems, and allows us to treat in a unified way a variety 
of situations, including systems with symmetry, morphisms, reduction, and 
nonlinearly constrained systems. Various examples illustrate the results. 


Contents 

1. Introduction 

2. Preliminaries 

2.1. Lie algebroids 

2.2. Exterior differential 

2.3. Morphisms 

2.4. Prolongation of a fibered manifold with respect to a Lie algebroid 

2.5. Prolongation of a map 

2.6. Lagrangian Mechanics 

3. Linearly constrained Lagrangian systems 

3.1. Lagrange-d’Alembert equations in local coordinates 

3.2. Solution of Lagrange-d’Alembert equations 

3.3. Projectors 

3.4. The distributional approach 

3.5. The nonholonomic bracket 

4. Morphisms and reduction 

4.1. Reduction of the free dynamics 

4.2. Reduction of the constrained dynamics 

4.3. Reduction by stages 

5. The momentum equation 

5.1. Unconstrained case 

5.2. Constrained case 

6. Examples 

6.1. Nonholonomic Lagrangian systems on Lie algebras 
The Suslov system 


2 

4 

4 

5 

5 

6 

7 

8 
I 

14 

15 

17 

18 
18 

19 

20 
21 

23 

25 

25 

26 
28 
28 
25 


Date : April 30, 2008. 

2000 Mathematics Subject Classification. 70F25, 70H03, 70H33, 37J60, 53D17. 

Key words and phrases. Nonholonomic Mechanics, Lagrange-d’Alembert equations, Lie alge¬ 
broids, symmetry, reduction. 

This work has been partially supported by Spanish Ministry of Education and Culture grants 
MTM2004-7832, BFM2003-01319, MTM2006-03322 and BFM2003-02532. J. Cortes was partially 
supported by faculty research funds granted by the University of California, Santa Cruz. 


1 







2 


J. CORTES, M. DE LEON, .7. C. MARRERO, AND E. MARTINEZ 


The Chaplygin sleigh 

6.2. Nonholonomic LR systems and right action Lie algebroids 
The Veselova system 

6.3. Semidirect product symmetry and left action Lie algebroids 

6.4. Chaplygin-type systems 
Chaplygin systems and Atiyah algebroids 

7. Nonlinearly constrained Lagrangian systems 

7.1. Dynamics in local coordinates 

7.2. Lagrange-d’Alembert solutions and nonholonomic bracket 

7.3. Morphisms and reduction 

7.4. Example: a ball rolling on a rotating table 

8. Conclusions and outlook 
References 


30 

31 
33 
35 
39 
44 

43 

44 

45 
47 
49 
53 
53 


1. Introduction 

The category of Lie algebroids has proved useful to formulate problems in ap¬ 
plied mathematics, algebraic topology, and differential geometry. In the context 
of Mechanics, an ambitious program was proposed in [SS] in order to develop for¬ 
mulations of the dynamical behavior of Lagrangian and Hamiltonian systems on 
Lie algebroids and discrete mechanics on Lie groupoids. In the last years, this 
program has been actively developed by many authors, and as a result, a powerful 
mathematical structure is emerging. 

The main feature of the Lie algebroid framework is its inclusive nature. Under 
the same umbrella, one can consider such disparate situations as systems with 
symmetry, systems evolving on semidirect products, Lagrangian and Hamiltonian 
systems on Lie algebras, and field theory equations (see 09j 02 for recent topical 
reviews illustrating this). The Lie algebroid approach to Mechanics builds on the 
particular structure of the tangent bundle to develop a geometric treatment of 
Lagrangian systems parallel to Klein’s formalism j22Jj26J. At the same time, the 
attention devoted to Lie algebroids from a purely geometrical viewpoint has led 
to an spectacular development of the field, e.g., see 0 mi m eqi and references 
therein. The merging of both perspectives has already provided mutual benefit, 
and will undoubtedly lead to important developments in the future. 

The other main theme of this paper are nonholonomic Lagrangian systems, i.e., 
systems subject to constraints involving the velocities. This topic is a classic sub¬ 
ject in Mathematics and Mechanics, dating back to the early times of Lagrange; 
a comprehensive list of classical references can be found in m- At the beginning 
of the nineties, the work E2 sparked a renewed interest in the geometric study 
of nonholonomic mechanical systems, with a special emphasis on symmetry as¬ 
pects. In the last years, several authors have extended the ideas and techniques of 
the geometrical treatment of unconstrained systems to the study of nonholonomic 
mechanical systems, see the recent monographs mm These include symplec- 
tic 021 snrni], Hamiltonian [53] , and Lagrangian approaches [221 SE] • the study of 
almost Poisson brackets 01121 [23, and symmetry and reduction of the dynam¬ 
ics [1010110107]. 

In this paper we develop a comprehensive treatment of nonholonomic systems 
on Lie algebroids. This class of systems was introduced in [20] when studying me¬ 
chanical control systems (see also 06] for a recent approach to mechanical systems 
on Lie algebroids subject to linear constraints). Here, we build on the geometry 
of Lie algebroids to identify suitable regularity conditions guaranteeing that the 
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nonholonomic system admits a unique solution. We develop a projection proce¬ 
dures to obtain the constrained dynamics as a modification of the unconstrained 
one, and define an almost-Poisson nonholonomic bracket. We show that many of 
the properties that standard nonholonomic systems enjoy have their counterpart 
in the proposed setup. As important examples, we highlight that the analysis here 
provides a natural interpretation for the use of pseudo-coordinates techniques and 
lends itself to the treatment of constrained systems with symmetry, following the 
ideas developed in [HI E3]- We carefully examine the reduction procedure for this 
class of systems, paying special attention to the evolution of the momentum map. 

From a methodological point of view, the approach taken in the paper has enor¬ 
mous advantages. This fact must mainly be attributed to the inclusive nature of Lie 
algebroids. Usually, the results on nonholonomic systems available in the literature 
are restricted to a particular class of nonholonomic systems, or to a specific context. 
However, as illustrated in Table [lj many different nonholonomic systems fit under 
the Lie algebroid framework, and this has the important consequence of making the 
results proved here widely applicable. With the aim of illustrating this breadth, 
we consider various examples throughout the paper, including the Suslov problem, 
the Chaplygin sleigh, the Veselova system, Chaplygin Gyro-type systems, the two¬ 
wheeled planar mobile robot, and a ball rolling on a rotating table. We envision 
that future developments within the proposed framework will have a broad impact 
in nonholonomic mechanics. In the course of the preparation of this manuscript, 
the recent research efforts mm were brought to our attention. These references, 
similar in spirit to the present work, deal with nonholonomic Lagrangian systems 
and focus on the reduction of Lie algebroid structures under symmetry. 


Nonholonomic 
Lagrangian system 

Lie algebroid 

Dynamics 

Example 

Standard 

Tangent bundle 

Lagrande-d’Alembert 

Rolling disk [57] 

On a Lie algebra 

Lie algebra 

Euler-Poincare-Suslov 

Chaplygin 
sleigh fTS] 

Nonholonomic LR 
systems 

Right action 

Lie algebroid 

Reduced Poincare- 

Chetaev 

Veselova problem 

m 

Nonholonomic systems 
with semidirect prod¬ 
uct symmetry 

Left action 

Lie algebroid 

Nonholonomic Euler- 
Poincare with an ad- 
vected parameter 

Chaplygin’s gyro 

Symmetry-invariant 

Atiyah algebroid 

Nonholonomic 

Lagrange-Poincare 

Snakeboard 


Table 1. The Lie algebroid framework embraces different classes 
of nonholonomic systems. 


The paper is organized as follows. In Section [2] we collect some preliminary no¬ 
tions and geometric objects on Lie algebroids, including differential calculus, mor- 
phisms and prolongations. We also describe classical Lagrangian systems within the 
formalism of Lie algebroids. In Section 02 we introduce the class of nonholonomic 
Lagrangian systems subject to linear constraints, given by a regular Lagrangian 
L : E —> R on the Lie algebroid r : E —> M and a constraint subbundle D 
of E. We show that the known results in Mechanics for these systems also hold 
in the context of Lie algebroids. In particular, drawing analogies with d’Alembert 
principle, we derive the Lagrange-d’Alembert equations of motion, prove the con¬ 
servation of energy and state a Noether’s theorem. We also derive local expressions 
for the dynamics of nonholonomic Lagrangian systems, which are further simplified 
by the choice of a convenient basis of D. As an illustration, we consider the class of 
nonholonomic mechanical systems. For such systems, the Lagrangian L is the polar 
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form of a bundle metric on E minus a potential function on M. In Section m 
we perform the analysis of the existence and uniqueness of solutions of constrained 
systems on general Lie algebroids, and extend the results in [a ® nni na no] for 
constrained systems evolving on tangent bundles. We obtain several characteriza¬ 
tions for the regularity of a nonholonomic system, and prove that a nonholonomic 
system of mechanical type is always regular. The constrained dynamics can be 
obtained by projecting the unconstrained dynamics in two different ways. Under 
the first projection, we develop a distributional approach analogous to that in [2], 
see also [56]. Using the second projection, we introduce the nonholonomic bracket. 
The evolution of any observable can be measured by computing its bracket with 
the energy of the system. Section B] is devoted to studying the reduction of the 
dynamics under symmetry. Our approach follows the ideas developed in m , who 
defined a minimal subcategory of the category of Lie algebroids which is stable un¬ 
der Lagrangian reduction. We study the behavior of the different geometric objects 
introduced under morphisms of Lie algebroids, and show that fiberwise surjective 
morphisms induce consistent reductions of the dynamics. This result covers, but 
does not reduce to, the usual case of reduction of the dynamics by a symmetry 
group. In accordance with the philosophy of the paper, we study first the uncon¬ 
strained dynamics case, and obtain later the results for the constrained dynamics 
using projections. A (Poisson) reduction by stages procedure can also be developed 
within this formalism. It should be noticed that the reduction under the presence 
of a Lie group of symmetries G is performed in two steps: first we reduce by a 
normal subgroup N of G, and then by the residual group. In Section [5] we prove 
a general version of the momentum equation introduced in [1] . In Section El we 
show some interesting examples and in Section 0 we extend some of the results 
previously obtained for linear constraints to the case of nonlinear constraints. The 
paper ends with our conclusions and a description of future research directions. 

2. Preliminaries 

In this section we recall some well-known facts concerning the geometry of Lie 
algebroids. We refer the reader to 01113 03 for details about Lie groupoids, Lie 
algebroids and their role in differential geometry. 

2.1. Lie algebroids. Let M be an n-dimensional manifold and let t : E — > M be 
a vector bundle. A vector bundle map p: E —> TM over the identity is called an 
anchor map. The vector bundle E together with an anchor map p is said to be 
an anchored vector bundle (see my a structure of Lie algebroid on E is 
given by a Lie algebra structure on the C°°(M)-module of sections of the bundle, 
(Sec(T), together with an anchor map, satisfying the compatibility condition 

WJv] = fW,v] + {p{<?)f)v- 

Here / is a smooth function on M, cr, p are sections of E and p(a) denotes the vector 
field on M given by p(a)(m) = p(a(m)). From the compatibility condition and the 
Jacobi identity, it follows that the map cr > p[cr) is a Lie algebra homomorphism 
from the set of sections of E , Sec(-E'), to the set of vector fields on M, X(M). 

In what concerns Mechanics, it is convenient to think of a Lie algebroid p: E — v 
TM, and more generally an anchored vector bundle, as a substitute of the tangent 
bundle of M. In this way, one regards an element a of £ as a generalized velocity, 
and the actual velocity v is obtained when applying the anchor to a, i.e., v = p(a). 
A curve a: [fo,fi] —► E is said to be admissible if m(t) = p(a(t)), where m(t) = 
T(a(t)) is the base curve. We will denote by Adm(T) the space of admissible curves 
on E. 
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Given local coordinates ( x 1 ) in the base manifold M and a local basis {e a } of 
sections of E , we have local coordinates ( x l ,y a ) in if. If a £ E is an element in 
the fiber over m £ M, then we can write a = y a e a (m ) and thus the coordinates of 
a are ( m l ,y a ), where m l are the coordinates of the point m. The anchor map is 
locally determined by the local functions p l a on M defined by p(e a ) = p l a (d/dx l ). 
In addition, for a Lie algebroid, the Lie bracket is determined by the functions 
defined by [e a ,ep] = C^^ej. The functions p l a and C(L are called the structure 
functions of the Lie algebroid in this coordinate system. They satisfy the following 
relations 


P 


o 

Ot 


d Pp 

dxi 





and 


cyclic(a:,/3, 7 ) 


dC ; 


01 


dx i 


rif-t 


= 0 , 


which are called the structure equations of the Lie algebroid. 


2.2. Exterior differential. The anchor p allows to define the differential of a 
function on the base manifold with respect to an element a £ E. It is given by 

df(a) = p{a)f. 

It follows that the differential of / at the point m £ M is an element of E 
Moreover, a structure of Lie algebroid on E allows to extend the differential to 
sections of the bundle /\ p E, which will be called p-sections or just p-forms. If 
w £ Sec (/\ p E), then du £ Sec(/\ P+1 E) is defined by 

duj(o 0 , or,.. •, ct p ) = 5>l)V( CTi )(u;(ao, a p )) 

i 

i<j 

It follows that d is a cohomology operator, that is, d 2 = 0. Locally the exterior 
differential is determined by 

dx l = p z a e a and de 1 = — ^ e“ A . 

Throughout this paper, the symbol d will refer to the exterior differential on the Lie 
algebroid E and not to the ordinary exterior differential on a manifold. Of course, 
if E = TM, then both exterior differentials coincide. 

The usual Cartan calculus extends to the case of Lie algebroids (see j3HJ [55j). 
For every section a oi E we have a derivation i a (contraction) of degree —1 and a 
derivation d a = i a o d + d o i a (Lie derivative) of degree 0. Since d 2 = 0, we have 
that d a o d = d o d G . 


2.3. Morphisms. Let r: E —> M and t' : E' —> M' be two anchored vector bun¬ 
dles, with anchor maps p: E —+ TM and p': E' —*■ TM'. A vector bundle map 
<1>: E —> E' over a map ip: M —> M' is said to be admissible if it maps admis¬ 
sible curves onto admissible curves, or equivalently T<pop = p'o<&. If E and E' 
are Lie algebroids, then we say that $ is a morphism if <&*d6 = d<&*9 for every 
9 £ Sec(/\ E'). It is easy to see that morphisms are admissible maps. 

In the above expression, the pullback <f>*/3 of a p-form fd is defined by 

($*/3) m (ai, 02 ,..., a p ) = ($(ai), $(a 2 ),..., $(a p )), 

for every ai,..., a p £ E m . For a function / £ C°°(M') (i.e., for p = 0), we just set 
$*/ = f °V- 

Let ( x l ) and (x n ) be local coordinate systems on M and M\ respectively. Let 
{ e a } and {e' a } be local basis of sections of E and E\ respectively, and {e Q } and 
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{e'“} the corresponding dual basis. The bundle map $ is determined by the rela¬ 
tions = <j> l (x) and <f>*e ,a = 4>^e^ for certain local functions <j> 1 and cj>p on M. 

Then, <1> is admissible if and only if 


dtp 

dxi 


P'Wo 


The map <f> is a morphism of Lie algebroids if and only if, in addition to the 
admissibility condition above, one has 




' dx i 


Ps 


9<t% 

dx i 



In these expressions, p l a , C? are the local structure functions on E and p H a , C'T 
are the local structure functions on E'. 


2.4. Prolongation of a fibered manifold with respect to a Lie algebroid. 

Let 7r: P —> M be a fibered manifold with base manifold M. Thinking of E 
as a substitute of the tangent bundle of M, the tangent bundle of P is not the 
appropriate space to describe dynamical systems on P. This is clear if we note that 
the projection to M of a vector tangent to P is a vector tangent to M, and what 
one would like instead is an element of E, the ‘new’ tangent bundle of M. 

A space which takes into account this restriction is the E-tangent bundle of 
P, also called the prolongation of P with respect to E, which we denote by T E P 
(see [321 ED EU1533) - It is defined as the vector bundle t e : T E P — ■> P whose fiber 
at a point p £ P m is the vector space 


T E P = { (6, v) £ E m x T p P | p(b) = T p n{v) }. 


We will frequently use the redundant notation (p, b, v ) to denote the element (b, v) £ 
T e P. In this way, the map t e is just the projection onto the first factor. The 
anchor of T E P is the projection onto the third factor, that is, the map p 1 : T E P —> 
TP given by p [ (p, b, v) = v. The projection onto the second factor will be denoted 
by Tn: T E P —► E, and it is a vector bundle map over n. Explicitly T7r(p, b, v) = b. 

An element z £ T E P is said to be vertical if it projects to zero, that is Ttt(z) = 0. 
Therefore it is of the form (j>, 0, v), with v a vertical vector tangent to P at p. 

Given local coordinates (x 1 , u A ) on P and a local basis {e a } of sections of E, we 
can define a local basis {X ai V^} of sections of T E P by 


( 9 

X a {p) = [p,e a {Tr{p)),p l a — 


3 and 


If z = (p,b,v) is an element of T E P , with b = z a e a , then v is of the form v = 
P l a z °‘'^P + and we can write 

2 = z a X a (p) + v a Va(p)- 


Vertical elements are linear combinations of {VU}- 

The anchor map p 1 applied to a section Z of T E P with local expression Z = 
Z a X a + V a Va is the vector field on P whose coordinate expression is 


P\Z) 


,Z' 


— 

dx l 




d 

du A ' 


If E carries a Lie algebroid structure, then so does T E P. The associated Lie 
bracket can be easily defined in terms of projectable sections, so that Tn is a 
morphism of Lie algebroids. A section Z of T E P is said to be projectable if there 
exists a section o of E such that Tn o Z = <j o n. Equivalently, a section Z is 
projectable if and only if it is of the form Z(p) = (p, o(n(p)), X(p)) : for some 
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section a of E and some vector field X on E (which projects to p(<r)). The Lie 
bracket of two projectable sections Z\ and Zi is then given by 

[Zi,Z 2 \(p) = {p,[a 1 ,a 2 ]{m),[X 1 ,X 2 \(p)), p G P, m = n(p). 

It is easy to see that [Z\, Z^lp) is an element of T E P for every p £ P. Since 
any section of T E P can be locally written as a linear combination of projectable 
sections, the definition of the Lie bracket for arbitrary sections of T E P follows. 
The Lie brackets of the elements of the basis are 

[X cn Xp\=CT 0 X 1 , K,Vb]= 0 and [V A ,V B } = 0, 
and the exterior differential is determined by 

dx^p^X 01 , du A = V A , 

dX 1 = -^Cl 0 X a A X 0 , dV A = 0, 

where {X a , V A } is the dual basis corresponding to {X a , Va}- 


2.5. Prolongation of a map. Let dt: P —> P' be a fibered map from the hbered 
manifold 7r: P — * M to the hbered manifold n': P' —► M' over a map p: M —* 
M’. Let $: E —> E' be an admissible map from r: E — ► M to t': E' —► M' 
over the same map p. The prolongation of d> with respect to T is the mapping 
T *'$ : T E P -> T E 'P' dehned by 

r*^{pAv) = (y{p),m,(T^){v)). 

It is clear from the definition that is a vector bundle map from t e : T E P —> P 
to t e , : T e P' —> P' over H/. Moreover, in [53] it is proved the following result. 

Proposition 2.1. The map T $ ’F is an admissible map. Moreover, Tis a 
morphism of Lie algebroids if and only if $ is a morphism of Lie algebroids. 


Given local coordinate systems (ad) on M and ( x ' 1 ) on M ', local adapted co¬ 
ordinates (x 1 ,u A ) on P and (x ,l ,u ,A ) on P' and a local basis of sections {e a } 
of E and {e^} of E ', the maps and T are determined by d>*e ,a = Q'Ze 0 and 
ty(x,u) = ((f> l (x), ip A (x, u)). Then the action of T $ \1/ is given by 


yr^^yx' 


= 


yr 9 ^yv ,A = p] 


dip' 


■X° 


d 


■V h 


dx i " ' du B 

We hnally mention that the composition of prolongation maps is the prolongation 
of the composition. Indeed, let 'F' be another bundle map from 7r': P' —* M' to 
another bundle 7r" : P" —► M" and be another admissible map from t' : E' —► M' 
to t" : E" —> M" both over the same base map. Since <f> and are admissible 
maps then so is o $, and thus we can define the prolongation of o f with 
respect to <f>' o d>. We have that T 9 o4> (iP' o T) = ( T $ T') o (T^'F). 

In the particular case when the bundles P and P' are just P = E and P' = E ', 
whenever we have an admissible map $: E —► E' we can define the prolongation 
of <I> along $ itself, by T 9 ${a, b, v) = (d>(a), <&(&), T$>(v)). From the result above, 
we have that T $ >F is a Lie algebroid morphism if and only if $ is a Lie algebroid 
morphism. In coordinates we obtain 


(T 4, $)*T , “ = &pX 0 , 

(T $ d>)*V'“ = p^-^y^X 0 + <t>%V 0 , 


where (x*,i/ 7 ) are the corresponding fibred coordinates on E. From this expression 
it is clear that is fiberwise surjective if and only if <F is fiberwise surjective. 
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2.6. Lagrangian Mechanics. In m (see also [H]) a geometric formalism for 
Lagrangian Mechanics on Lie algebroids was defined. Such a formalism is similar 
to Klein’s formalism [36] in standard Lagrangian mechanics and it is developed in 
the prolongation T E E of a Lie algebroid E over itself. The canonical geometrical 
structures defined on T e E are the following: 

• The vertical lift £ v : t*E — > T e E given by £ v (a, b ) = (a, 0, b%), where b% 
is the vector tangent to the curve a + tb at t = 0, 

• The vertical endomorphism S: T e E —> T e E defined as follows: 

S(a, 6, v) = C(a, b) = (a, 0, b v a ), 

• The Liouville section which is the vertical section corresponding to the 
Liouville dilation vector field: 


A(a) = C(a,a) = (a, 0,a%). 

A section T of T e E is said to be a SODE section if ST = A. 

Given a Lagrangian function L £ C°°(E) we define the Cartan 1-form 0l and 
the Cartan 2- form as the forms on T e E given by 

0L = S*(dL) and lol = —dO^. (2-1) 

The real function El on E defined by El = d&L — L is the energy function of 
the Lagrangian system. 

By a solution of the Lagrangian system (a solution of the Euler-Lagrange 
equations) we mean a SODE section T of T e E such that 

ir<-L>L — AEl = 0 . ( 2 - 2 ) 


The local expressions for the vertical endomorphism, the Liouville section, the 
Cartan 2-form and the Lagrangian energy are 

(2.3) 

(2.4) 


SX a = Va, SVa = 0, for all a, 


A = y a v a , 


lu l = 




dy a dyP 


2 \dx l dy a 1 dx z dyP 


^ + (2.5) 

( 2 . 6 ) 


Thus, a SODE T is a section of the form 


r = y a x a + rv a . 


The SODE T is a solution of the Euler-Lagrange equations if and only if the functions 
/“ satisfy the linear equations 


d 2 L 

dy&dy c 


■r + 


d 2 L 

dx i dy c 


dL 


dL 


PW "Pa^ = °, a. 


dx i 


(2.7) 


The Euler-Lagrange differential equations are the differential equations for 
the integral curves of the vector field / 0 1 (T), where the section T is the solution of 
the Euler-Lagrange equations. Thus, these equations may be written as 


z = Pair, 


d dL t dL dL „ g 
df “ Pa dX* + d^ Ca ^ y 


= 0 . 


In other words, if 5L : Adm(E) —> E* is the Euler-Lagrange operator , which 
locally reads 


d , dL 


dL , dL 
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where {e“} is the dual basis of {e a }, then the Euler-Lagrange differential equations 
read 

SL = 0. 

The function L is said to be regular Lagrangian if u>l is regular at every point 
as a bilinear map. In such a case, there exists a unique section T^ of T E E which 
satisfies the equation 

ir h u)L — dE^ = 0 . 

Note that from (12.3[) . (12.41) . (|2.5D and (|2.6D . it follows that 

isx^L = ~S*(ixU)L), = —S*(dE L ), (2.8) 

for X £ Sec {T e E). Thus, using (12.81) . we deduce that 


isv L UL = i/\UL 


which implies that T^ is a SODE section. Therefore, for a regular Lagrangian func¬ 
tion L we will say that the dynamical equations (|2.2I) are just the Euler-Lagrange 
equations. 

On the other hand, the vertical distribution is isotropic with respect to wl, 
see |33] • This fact implies that the contraction of ujl with a vertical vector is a 
semibasic form. This property allows us to define a symmetric 2-tensor G L along 
r by 

G L a (b,c)=co L (b,c v a ), _ (2.9) 

where b is any element in T E E which projects to &, i.e., Tr(&) = &, and a £ E. In 
coordinates G L = W a pe a ® e@, where the matrix W a p is given by 


g 2 t 

Wa P = dy a dyE 

It is easy to see that the Lagrangian L is regular if and only the matrix W is regular 
at every point, that is, if the tensor G L is regular at every point. By the kernel of 
G L at a point a we mean the vector space 


ICer Ga = { b £ E r{a) \ (6, c) = 0 for all c e E r{a) }. 

In the case of a regular Lagrangian, the Cartan 2-section ujl is symplectic (non¬ 
degenerate and d-closed) and the vertical subbundle is Lagrangian. It follows that a 
1-form is semi-basic if and only if it is the contraction of ujl with a vertical element. 

Finally, we mention that the complete lift a c of a section a £ Sec(If) is the 
section of T e E characterized by the two following properties: 

(1) projects to a, i.e., Ttoo c =oot, 

(2) d a cjl — dcr/i, 

where by a £ C°°(E) we denote the linear function associated to a 1-section a £ 
Sec (E*). Note that 

T SODE section, a £ Sec(if) => ^[cr 0 , T] = 0, (2.11) 

57 = 0, ct G Sec(E) => 5[cr c ,7] = 0. (2.12) 


3. Linearly constrained Lagrangian systems 

Nonholonomic systems on Lie algebroids were introduced in [5D]. This class of 
systems includes, as particular cases, standard nonholonomic systems defined on 
the tangent bundle of a manifold and systems obtained by the reduction of the 
action of a symmetry group. The situation is similar to the non-constrained case, 
where the general equation 5L = 0 comprises as particular cases the standard 
Lagrangian Mechanics, Lagrangian Mechanics with symmetry, Lagrangian systems 
with holonomic constraints, systems on semi-direct products and systems evolving 
on Lie algebras, see e.g., m ■ 
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We start with a free Lagrangian system on a Lie algebroid E. As mentioned 
above, these two objects can describe a wide class of systems. Now, we plug in 
some nonholonomic linear constraints described by a subbundle D of the bundle 
E of admissible directions. If we impose to the solution curves a(t ) the condition 
to stay on the manifold D , we arrive at the equations <5L a ( t ) = A (t) and aft) G D, 
where the constraint force A (t) G E*r a (t)) f° determined. In the tangent bundle 
geometry case ( E = TM), the d’Alembert principle establishes that the mechanical 
work done by the constraint forces vanishes, which implies that A takes values in 
the annihilator of the constraint manifold D. Therefore, in the case of a general 
Lie algebroid, the natural equations one should pose are (see EUj ) 

^a(t) 6 ^r(a(t)) an d £ D. 

In more explicit terms, we look for curves aft) on E such that 

- they are admissible, p(aft)) = rhft ), where m = t o a, 

- they stay in D, aft) G D m ( t ), 

- there exists A (t) G D^tt) such that <5L 0 ( t ) = Xft). 

If aft) is one of such curves, then (a(t),a(t)) is a curve in T e E. Moreover, since 
aft) is in U, we have aft) is tangent to D, that is, (aft), aft)) G T d D. Under 
some regularity conditions (to be made precise later on), we may assume that 
the above curves are integral curves of a section T, which as a consequence will 
be a SODE section taking values in T D D. Based on these arguments, we may 
reformulate geometrically our problem as the search for a SODE T (defined at least 

on a neighborhood of D) satisfying (ir^L — dE^a G D°( a \ and T(a) G T E D, at 
every point a G D. In the above expression D° is the pullback of D° to T e E , that 
is, a G D^ia.) ^ an d only if there exists A G D°^ such that a = A o T a r. 

Definition 3.1. A nonholonomically constrained Lagrangian system on a Lie 

algebroid E is a pair (L, D), where L is a smooth function on E, the Lagrangian , 
and i: D E is a smooth subbundle of E, known as the constraint subbundle. 
By a solution of the nonholonomically constrained Lagrangian system (L, D) we 
mean a SODE section T G T e E which satisfies the Lagrange-d’Alembert equa¬ 
tions _ 

(iraJL — dE L )\ D G Sec(£>°), 

T| D G Sec (T d D). 

With a slight abuse of language, we will interchangeably refer to a solution of the 
constrained Lagrangian system as a section or the collection of its corresponding 
integral curves. The restriction of the projection r: E —>■ M to D will be denoted 
by 7T , that is, 7r = t\d '■ D —> M. 

Remark 3.2 (Domain of definition of solutions of the Lagrange-d’Alembert equa¬ 
tions). We want to stress that a solution of the Lagrange-d’Alembert equations 
needs to be defined only over D 1 but for practical purposes we consider it extended 
to E (or just to a neighborhood of D in E). We will not make any notational dis¬ 
tinction between a solution on D and any of its extensions. Solutions which coincide 
on D will be considered as equal. See [301 03 for a more in-depth discussion. In 
accordance with this convention, by a SODE on D we mean a section of T D D which 
is the restriction to D of some SODE defined in a neighborhood of D. Alternatively, 
a SODE on D is a section T of T d D such that Tr(T(a)) = a for every a G D. • 

Remark 3.3 (Holonomic constraints). A nonholonomically constrained Lagrangian 
system (L, D) on a Lie algebroid E is said to be holonomic if D is a Lie subalge- 
broid of E. This means that [X, Y] G Sec(D), for X, Y G Sec(D). Thus, the real 
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function Ld = L\ D : D —> M defines an unconstrained (free) Lagrangian system 
on the Lie algebroid D. Moreover, it is easy to prove that L o Ad = A o i and 
1 o Sd = Sol, where T = T l i : T d D —> T e E is the prolongation of the Lie 
algebroid morphism i: D E and Ad (respectively, Sd) is the Liouville section 
(respectively, the vertical endomorphism) of the Lie algebroid T d D. Therefore, 
since L o i = Ld , we deduce that 

1*(0 l ) = 0 Ld , 1*{u l )=u L d , T{dE L )=dE LD . 

Consequently, if T is a SODE section of T e E , a,b G D, ( b , X ) € T e E and (6, Y) G 
T e D then 

(twl - dE L )(a)(b,X) = {ir lD UL D ~ dE LD )(a)(b,Y) + {ivto L - dE L )(a)(Q,Z), 

(0, Z) being a vertical element of T e E. 

Now, using (12.81) . we have that (ir^L — dEi,){a){ 0, Z) = 0 which implies that 

(irUL - dE L )(a)(b, X) = (ir lD v LD - dE Lo )(a)(b,Y). 

The above facts prove that a SODE section T of T e E is a solution of the holonomic 
Lagrangian system (L, D) on E if and only if T|d is a solution of the Euler-Lagrange 
equations for the (unconstrained) Lagrangian function Ld on the Lie algebroid D. 
In other words, the holonomic Lagrangian system ( L,D) on E may be considered 
as an unconstrained (free) Lagrangian system on the Lie algebroid D. • 

Next, suppose that (L, D ) is a nonholonomically constrained Lagrangian system 
on the Lie algebroid E. Then, the different spaces wc will consider are shown in 
the following commutative diagram 


TM TM 



M M 

As an intermediate space in our analysis of the regularity of the constrained 
systems, we will also consider T E D , the .E-tangent to D. The main difference 
between T E D and T D D is that the former has a natural Lie algebroid structure 
while the later does not. 

The following two results are immediate consequences of the above form of the 
Lagrange-d’Alembert equations. 

Theorem 3.4 (Conservation of energy). If (L,D) is a constrained Lagrangian 
system and T is a solution of the dynamics, then dr El =0 (on D). 

Proof. Indeed, for every a £ D, we have T(a) € T E D, so that Tr(T(a)) 6 D. 
Therefore irD° = 0 and contracting 0 = ir{ir^>L — dEjf) = —drEL at every point 
in D. □ 
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Theorem 3.5 (Noether’s theorem). Let ( L,D) be a constrained Lagrangian system 
which admits a unique SODE F solution of the dynamics. If a is a section of D 
such that there exists a function f £ C°°(M) satisfying 

d a cL = /, 

then the function F = (6r,cr c ) — f is a constant of the motion, that is, drF = 0 
(on D). 

Proof. Using that 9l{ r) = d&(L), we obtain i a c (ira>L — dE^) = i a c (— drOr+dL) = 
d a cL — dr ( 0L: & c ) + 9l[L, a c ] and, since [r, <r c ] is vertical, we deduce 

i a c(i r ujL — dEjf) = d a cL — dr ( 9l , <x c ) ■ 

Thus, taking into account that i a cD° = 0, we get 0 = dr{(0L,o' c ) — f) = —drF. □ 

Example 3.6 (Mechanical systems with nonholonomic constraints). Let Q : E Xjh 
E —► R be a bundle metric on E. The Levi- Civita connection V s is determined 
by the formula 

2^(V^?7,C) =p(<r)(&(r}, 0) +p(ii)(G{<J,Q) - P(0(G(v,cr)) 

+ G(e, [C,V])+G{V, [C,cr]) - G((, [t7 , cr]), 
for cr, 77 , £ £ Sec (E). The coefficients of the connection V e are given by 

r /? 7 = \G av {W,fci\ + [ft,P\ + \P,rA), 

where Q av are the coefficients of the metric Q, (G av ) is the inverse matrix of (Gau) 
and 

Using the covariant derivative induced by V s , one may introduce the notion of a 
geodesic of as follows. An admissible curve a : I —> E is said to be a geodesic if 
V® t ,a(t) = 0, for all t £ I. In local coordinates, the conditions for being a geodesic 
read 

+ ^( r «/3 + r iL )«“ a/3 = °> for a11 7 - 

The geodesics are the integral curves of a SODE section V v c of T e E, which is locally 
given by 

r vs = y UT 7 -i(r^ + ry y Vv 7 . 

r vS is called the geodesic flow (for more details, see m)- 

The class of systems that were considered in detail in [2D] is that of mechanical 
systems with nonholonomic constraint ;f|. The Lagrangian function L is of 
mechanical type, i.e., it is of the form 

L(a) = ^G(a, a) - V (r(a)), a £ E, 

with V a function on M. 

The Euler-Lagrange section for the unconstrained system can be written as 

Fl = r vs - (grad 6 V) v . 

In this expression, by grad e V we mean the section of E such that (dV(m),a) = 
C?(grad e V(m),a), for all m £ M and all a £ E m , and where we remind that d is 


^In fact, in m, we considered controlled mechanical systems with nonholonomic constraints, 
that is, mechanical systems evolving on Lie algebroids and subject to some external control forces. 
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the differential in the Lie algebroid. The Euler-Lagrange differential equations can 
be written as 

& = piv a , 

1 r)V (3-2) 

Alternatively, one can describe the dynamical behavior of the mechanical control 
system by means of an equation on E via the covariant derivative. An admissible 
curve a : 1 i—» a(t) with base curve t > m(t) is a solution of the system (13.21) if and 
only if 

Vf (t) a(£) + grad g V(m(t)) = 0. (3.3) 

Note that 

0MO)(Vf (t) a(f) + grad e V(m(t)), b) = SL(a(t))(b ), for b £ E m{t) . 

If this mechanical control system is subject to the constraints determined by a 
subbundle D of E, we can do the following. Consider the orthogonal decomposition 
E = D © D- 1 , and the associated orthogonal projectors P : E —► D, Q : E —> D A . 
Using the fact that Q(P-,-) = one can write the Lagrange-d’Alembert 

equations in the form 

P ( V f(t) a W) + -P(gradg V{m(t))) = 0, Q(a) = 0. 

A specially well-suited form of these equations makes use of the constrained con¬ 
nection V defined by \7 a r] = P(V®?y) + V^iQrj). In terms of V, we can rewrite 
this equation as V a ( t )a(£) + P(grad e V(m(t))) = 0, Q{a) = 0, where we have used 
the fact that the connection V restricts to the subbundle D. 

Moreover, following the ideas in [J5|, we proved in (2U] that the subbundle D 
is geodesically invariant for the connection V, that is, any integral curve of the 
spray Ty associated with V starting from a point in D is entirely contained in D. 
Since the terms coming from the potential V also belongs to D, we have that the 
constrained equations of motion can be simply stated as 

V a (t)Ci(t) + P{ grad g = °> a (°) e D - ( 3 - 4 ) 

Note that one can write the constrained equations of the motion as follows 
a{t) = p 1 (Ty (a(i)) - P(grad a V) v {a{t))) 

and that the restriction to D of the vector field / 0 1 (Ty — P(grad e V")") is tangent 
to D. 

The coordinate expression of equations (ECU) is greatly simplified if we take a 
basis {e a } = {e a , e^} of E adapted to the orthogonal decomposition E = D ® £>-*-, 
i.e., D = span{e a }, PU = spanje^}. Denoting by (y a ) = {y a ,y A ) the induced 
coordinates, the constraint equations Q(a) = 0 just read y A = 0. The differential 
equations of the motion are then 

X 1 = p l a y a , 

y a = -\{T a hc + t a cb ) yY-G ab pi^-, 

y A = 0 , 

where are the connection coefficients of the constrained connection V. • 

In the above example the dynamics exists and is completely determined whatever 
the (linear) constraints are. As we will see in Section UlZI this property is lost in 
the general case. 
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3.1. Lagrange-d’Alembert equations in local coordinates. Let us analyze 
the form of the Lagrange-d’Alembert equations in local coordinates. Following the 
example above, let us choose a special coordinate system adapted to the structure 
of the problem as follows. We consider local coordinates (x l ) on an open set U 
of M and we take a basis {e a } of local sections of D and complete it to a basis 
{e a , e^} of local sections of E (both defined on the open U). In this way, we have 
coordinates (ad, y a , y A ) on E. In this set of coordinates, the constraints imposed by 
the submanifold D C E are simply y A = 0. If {e a , e A } is the dual basis of {e Q , e^}, 
then a basis for the annihilator D° of D is {e* 4 } and a basis for D° is X A . 

An element z of T E D is of the form 2 = u a X a + z a V a = u a X a + u a Xa + z a V a , 
that is, the component Va vanishes since p 1 (z) is a vector tangent to the manifold 
D with equations y A = 0. The projection of z to E is Tt(z) = u a e a + u A eA , so 
that the element 2 is in T D D if and only if u A = 0. In other words, an element in 
T d D is of the form 2 = u a X a + z a V a - 

Let us find the local expression of the Lagrange-d’Alembert equations in these 
coordinates. We consider a section T such that G Sec (T d D), which is therefore 
of the form T = g a X a + f a V a ■ From the local expression (|2.5p of the Cartan 2-form 
and the local expression (12.tip of the energy function, we get 


0 = (* r wL - dE L ,V a ) = -y B 


d 2 L 

dy a dy B 


(: y b 


9 b ) 


d 2 L 

dy a dy b ' 


If we assume that the Lagrangian L is regular, when we evaluate at y A = 0, we 
have that g a = y a and thus T is a SODE. Moreover, contracting with X a , after a 
few calculations we get 


0 = {itUl - dE L ,X a ) 



^Lra n.0 _ J d-k 
dy'y a0y Pa dx i 


so that (again after evaluation at y A = 0 ), the functions f a are solution of the 
linear equations 


d 2 L 

dy b dy a 


f b + 


d 2 L 

dx l dy a 


PbV 


v b - P l — = 0 
d y i abV Pa dx i 


(3.5) 


where all the partial derivatives of the Lagrangian are to be evaluated on y A = 0. 

As a consequence, we get that there exists a unique solution of the Lagrange- 
d’Alembert equations if and only if the matrix 


c ab {x i ,y c ) = -^— s {x\y c , 0) (3.6) 

is regular. Notice that C a b is a submatrix of W a p, evaluated at y A = 0 and that, as 
we know, if L is of mechanical type then the Lagrange-d’Alembert equations have 
a unique solution. The differential equations for the integral curves of the vector 
field p 4 (r) are the Lagrange-d’Alembert differential equations, which read 


& = , 

±(*L\ dL_0u_ = 

dt\dy a ) dy~< ab!J Pa dx i 

y A = 0. 


(3.7) 


Finally, notice that the contraction with Xa just gives the components Xa = 
(irbOL — dEL, Xa) of the constraint forces A = A a& A - 
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Remark 3.7 (Equations in terms of the constrained Lagrangian). In some occa¬ 
sions, it is useful to write the equations in the form 


= PaV a , 

d_ f 0L\ SL b _ idL = 8L A 
dt \dy a ) + dy^ abV Pa 8x i dy ACabV 

y A = 0 , 


(3.8) 


where, on the left-hand side of the second equation, all the derivatives can be 
calculated from the value of the Lagrangian on the constraint submanifold D. In 
other words, we can substitute L by the constrained Lagrangian L c defined by 
L c {x\y a ) = L{x\y\ 0). 


Remark 3.8 (Lagrange-d’Alembert equations in quasicoordinates). A particular 
case of this construction is given by constrained systems defined in the standard 
Lie algebroid tm'- TM —> M. In this case, the equations (1X71) are the Lagrange- 
d’Alembert equations written in quasicoordinates, where C are the so-called 
Hamel’s transpositional symbols, which obviously are nothing but the structure 
coefficients (in the Cartan’s sense) of the moving frame { e a }, see e.g., [2a EU. • 


3.2. Solution of Lagrange-d’Alembert equations. 


Assumption 3.9. In what follows, we will assume that the Lagrangian L is regular 
at least in a neighborhood of D. 


Let us now perform a precise global analysis of the existence and uniqueness of 
the solution of Lagrange-d’Alembert equations. 

Definition 3.10. A constrained Lagrangian system ( L,D ) is said to be regular if 
the Lagrange-d’Alembert equations have a unique solution. 


In order to characterize geometrically those nonholonomic systems which are 
regular, we define the tensor G L,D as the restriction of G L to D, that is, G^’ D (b, c) = 
G^(b,c) for every a £ D and every b, c £ D r r a y In coordinates adapted to D, we 
have that the local expression of G L ’ D is G L,D = C a be a <8> e b where the matrix C a b is 
given by equation 

A second important geometric object is the subbundle F C T e E\d D whose 
fiber at the point a £ D is F a = uf 1 (D°^). More explicitly, 

F a = {z £ T e E | exists ( £ D°^ s.t. lol{z,u ) = ((,Tt(u)} for all u £ T e E}. 

From the definition, it is clear that the rank of F is rank(F) = rank(D°) = 
rank(A') — rank(D). 

Finally, we also consider the subbundle ( T d D) ± C T e E\d — >■ D, the orthog¬ 
onal to T d D with respect to the symplectic form u>l- The rank of (T d D) ± is 
rank {T D D) 1 - = rank(T £ '£’) — rank(T D D) = 2(rank(£() — rank(D)) = 2rank(D°). 

The relation among these three objects is described by the following result. 


Lemma 3.11. The following properties are satisfied: 

(1) The elements in F are vertical. An element £ v (a,b) £ F n if and only if 
Ga(b, c) = 0 for allc£D T(a) . 

(2) (rV-Lnver(T*B) = ; 


Proof. (1) The elements in F are vertical because the elements in D° are semi-basic. 
If f v (a,b) £ F a then there exists £ G £>°. ^ such that wl(£ v (<z, b), u) = (( ,Tt(u )) 
for all u £ T e E. In terms of G L and writing c = Tt(u), the above equation 
reads —G^(b,c) = ((,c). By taking u £ Tt~ 1 (D), then c is in D and therefore 
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G L a (b,c) = 0 for all c G D T ^ a y Conversely, if G^(b,c) = 0 for all c G -D r ( 0 )> then 
the 1-form £ = —G F {b, ) is in D°( a y Therefore wl(£ v (o, b), u) = — G F {b, Tt(u)) = 
{(,Tr(n))), which is the condition for f v (a, 6) G F a . 

(2) The condition for a vertical element £ v (a, 6) to be in (T d D) ± is 6), w) = 

0 for all io G T a D D, or equivalently, G^(b 1 Tr(w)) = 0. The vector c = Tt(w) is an 
arbitrary element of -D r (a), so that the above condition reads G%(b, c) = 0, for all 
c G D r ( a y which is precisely the condition for f v (a, b) to be in F„. □ 

Theorem 3.12. The following properties are equivalent: 

(1) The constrained Lagrangian system ( L,D ) is regular, 

(2) Ker G L,D = {0}, 

(3) r E Dnf’ = {o}, 

(4) T D Dn(T D D) ± = {0}. 

Proof. [(l)o(2)] The equivalence between the first two conditions is clear from the 
local form of the Lagrange-d’Alembert equations (13.51) . since the coefficients of the 
unknowns /“ are precisely the components (13.61) of G L,D . 

[(2)4>(3)] (=>) Let a G D and consider an element z G T E D n F a . Since the 
elements of F are vertical, we have 2 = f v (a,b) for some b G E T ( a y Moreover, 
2 G T E D implies that b is an element in D T ^ a y On the other hand, if z = f v (a,b) 
is in F a , then Lemma 13.111 implies that G F (b,c) = 0 for all c G H r ( a )- Thus 
G ED (b, c) = 0 for all c G Z) T ( Q ), from where 6 = 0, and hence z = 0. 

(<=) Conversely, if for some a G D, there exists 6 G KerG ED with 6^0 then, 
using Lemma [3.Ill we deduce that z = f v {a, b) G F 0 Fj D n F a and z ^ 0. 

[(2)4>(4)] (=>) Let a G D and consider an element v G T E D 0 (T E D)^~, that 
is, ujl(v,w) = 0 for all w G T a ,J D. If we take w = £ v (a,b) for 6 G F > T ( a ) arbitrary, 
then we have ujl{v,^ v (a, 6)) = G ED (Tr(v),b) = 0 for all 6 G D T ( a y from where it 
follows that Tt(v) = 0. Thus v is vertical, v = £ v (a, c), for some c G D and then 
Wi({ v (a,c),w) = —G ed (c,Tt(w)) = 0 for all w G F FJ D. Therefore c = 0 and 
hence v = 0. 

(4=) Conversely, if for some a G D, there exists 6 G Ker G ED with 6^0, then 
0 / £ v ( a > b) G T E D fl {T E D) 1 -, because wl(£ v (o, 6), w) = G L ’ D (b,Tr(w)) = 0 for 
all w G T E D. □ 

In the case of a constrained mechanical system, the tensor G L is given by 
G F (b,c) = G T ( a ){b,c), so that it is positive definite at every point. Thus the re¬ 
striction to any subbundle D is also positive definite and hence regular. Thus, 
nonholonomic mechanical systems are always regular. 

Proposition 3.13. Conditions (3) and (4) in Theorem \3.12\ are equivalent, respec¬ 
tively, to 

(3’) T e E\ d = T e D®F, 

(4’) T e E\ d =T d D®(T d D) ± . 

Proof. The equivalence between (4) and (4’) is obvious, since we are assuming that 
the free Lagrangian is regular, i.e., ujl is symplectic. The equivalence of (3) and 
(3’) follows by computing the dimension of the corresponding spaces. The ranks of 
T e E , T e D and F are 

rank (T e E) = 2rank(G), 
rank(T' B H) = rank(£?) + rank(H), 
rank(P) = rank(D°) = rank(£?) — rank(D). 

Thus rank (T e E) = rank(T E D) + rank(F), and the result follows. □ 
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3.3. Projectors. We can express the constrained dynamical section in terms of 
the free dynamical section by projecting to the adequate space, either T E D or 
T D D , according to each of the above decompositions of T e E\b- Of course, both 
procedures give the same result. 

Projection to T E D. Assuming that the constrained system is regular, we have a 
direct sum decomposition 

T e E = T e D © F a . 

for every a £ D, where we recall that the subbundle F C T E D is defined by 
F = w£ 1 (D°), or equivalently D° = ujl{F). 

Let us denote by P and Q the complementary projectors defined by this decom¬ 
position, that is, 

P a : T e E -> T e D and Q a : T e E -► F a , for all a £ D. 

Then we have, 

Theorem 3.14. Let (L, D ) be a regular constrained Lagrangian system and let T^ 
be the solution of the free dynamics, i.e., ir L u>L = dE l- Then the solution of the 
constrained dynamics is the SODE T(l,d) obtained by projection = P(T l\d)- 

Proof. Indeed, if we write = Tl(o) — Q(Tl{o)) for a £ D, then we have 

*r(L,D)(») w i - dE L (a) = ir L ( a )^L ~ ^Q(r L (a))^L ~ dE L (a ) = -*Q( rr , (»))<*>£ G D ° T { a y 

which is an element of D°^ because Q(Tl(o)) is in F a . Moreover, since is a 
SODE and Q(T l) is vertical (since it is in F), we have that P(l,d) is also a SODE. □ 

We consider adapted local coordinates (x l , y a ,y A ) corresponding to the choice of 
an adapted basis of sections {e a ,eA}, where {e a } generate D. The annihilator D° 
of D is generated by {e A }, and thus D° is generated by {T' 4 }. A simple calculation 
shows that a basis { Za } of local sections of F is given by 

Za = Va - Q a A Va, (3.9) 

where Q a A = WAbC ab and C ab are the components of the inverse of the matrix C a b 
given by equation (EH). The local expression of the projector over F is then 

Q = Z a ®V a . 

If the expression of the free dynamical section T^ in this local coordinates is 

r L = y a x a + f a v a , 

(where /“ are given by equation (I2.7D ). then the expression of the constrained 
dynamical section is 

r (L,D) = y a Xa + (f a + f A Q a A )Va, 
where all the functions f a are evaluated at y A = 0. 

Projection to T D D. We have seen that the regularity condition for the constrained 
system (L, D) can be equivalently expressed by requiring that the subbundle T D D 
is a symplectic subbundle of (T E E,ojl). It follows that, for every a £ D, we have 
a direct sum decomposition 

T e E = T a D D © ( T a D D) E . 

Let us denote by P and Q the complementary projectors defined by this decompo¬ 
sition, that is, 

P a :T E E^T a D D and Q a : T E E - {T a ° , for all a £ D. 

Then, we have the following result: 
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Theorem 3.15. Let ( L,D ) be a regular constrained Lagrangian system and let 
be the solution of the free dynamics, i.e., ir L u>L = dE^. Then the solution of the 
constrained dynamics is the SODE T^m obtained by projection ^ = P(Tl\d). 

Proof. From Theorem 13.141 we have that the solution F^^q of the constrained 
dynamics is related to the free dynamics by T^Id = F + (3(r£|£>). Let us prove that 
Q{^l\d) takes values in (T D D) 1 - . Indeed, Q(Tl\d) takes values in F = (T D D) ± n 
Ver (T e E), so that, in particular, it takes values in ( T D D) J Thus, since F is a 
section of T D D, it follows that Tl\d = d) is a decomposition of Vl\d 

according to T e E\d = T D D ® (T d D) ± , which implies T(l,d) = P(Tl\d)- □ 

In adapted coordinates, a local basis of sections of ( T d D) ± is { Ya,Za }, where 
the sections Z A are given by and the sections Ya are 

Ya = -Fa - Q a A*a + C bc {M Ab - M ab Q a A )Vc, 

with M a p = ujL(X a , Xp). Therefore the expression of the projector onto (T D D) E 
is 

Q = Z A ®V A + Y a ®X a . 

Note that S(Y A ) = Z A - 

3.4. The distributional approach. The equations for the Lagrange-d’Alembert 
section T can be entirely written in terms of objects in the manifold T D D. Recall 
that T d D is not a Lie algebroid. In order to do this, define the 2-section io L ’ D 
as the restriction of wl to T d D. If ( L,D ) is regular, then T D D is a symplectic 
subbundle of (T E E,u>l)- From this, it follows that u> L,D is a symplectic section on 
that bundle. We also define e L,D to be the restriction of dE^ to T D D. Then, taking 
the restriction of Lagrange-d’Alembert equations to T D D, we get the following 
equation 

i r uj L ’ D = e L ’ D , (3.10) 

which uniquely determines the section F. Indeed, the unique solution T of the above 
equations is the solution of Lagrange-d’Alembert equations: if we denote by A the 
constraint force, we have for every u € T E D that 

uj L (T(a),u) - (dE L (a),u ) = (A(o),Tt(m)) = 0, 

where we have taken into account that Tt(u) € D and A(a) € D°. 

This approach, the so called distributional approach , was initiated by Bocha¬ 
rov and Vinogradov (see [67] ) and further developed by Sniatycki and coworkers j2[ 
1231E3■ Similar equations, within the framework of Lie algebroids, are the base of 
the theory proposed in ]56j . 

Remark 3.16 (Alternative description with T E D). One can also consider the re¬ 
striction to T E D, which is a Lie algebroid, but no further simplification is achieved 
by this. If Co is the restriction of u>l to T E D and e is the restriction of dE l to T E D, 
then the Lagrange-d’Alembert equations can be written in the form irio — E = A, 
where A is the restriction of the constraint force to T E D, which, in general, does 
not vanish. Also notice that the 2-form Co is closed but, in general, degenerated. • 

3.5. The nonholonomic bracket. Let f,g be two smooth functions on D and 
take arbitrary extensions to E denoted by the same letters (if there is no possibility 
of confusion). Suppose that Xf and X g are the Hamiltonian sections on T e E given 
respectively by 

ix f wl = df and i Xg u L =dg. 

We define the nonholonomic bracket of / and g as follows: 

{f,g}„ h =u L (P(X f ),P(X g )). 


(3.11) 
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Note that if /' is another extension of /, then (Xf — Xp)\ D is a section of ( T d D) ± 
and, thus, we deduce that (|3. 1 II) does not depend on the chosen extensions. The 
nonholonomic bracket is an almost-Poisson bracket, i.e., it is skew-symmetric, a 
derivation in each argument with respect to the usual product of functions and 
does not satisfy the Jacobi identity. 

In addition, one can prove the following formula 

f={f,E L } nh . ( 3 . 12 ) 

Indeed, we have 

/ = d r (L ,D)f = *r {L , D) d f = ir( L ' D) ix f UL 

= u L (x f ,r (LiD) ) = u L (x f ,P(r L )) 

= u L (P{x f ),p{r L )) = {f,E L } nh . 

Equation (13.121) implies once more the conservation of the energy (by the skew- 
symmetric character of the nonholonomic bracket). 

Alternatively, since T D D is an anchored vector bundle, one can take the function 
/ € C co (D) and its differential df £ Sec ((T d D)*). Since uj l,d is regular, we have a 
unique section Xf £ Sec (T D D) defined by ix f u> L,D = df. Then the nonholonomic 
bracket of two functions / and g is {/, g} n h = u L ’ D (Xf,X g ). Note that if / £ 
C°°(E) (resp. g £ C°°(E)) is an extension to E of / (resp., g ), then Xf = P{Xj) \r> 
(resp., X g = P(Xg) {D ). 

4. MORPHISMS AND REDUCTION 


One important advantage of dealing with Lagrangian systems evolving on Lie 
algebroids is that the reduction procedure can be naturally handled by considering 
morphisms of Lie algebroids, as it was already observed by Weinstein [55]. We 
study in this section the transformation laws of the different geometric objects in 
our theory and we apply these results to the study of the reduction theory. 

Proposition 4.1. Let 4>: E —> E’ be a morphism of Lie algebroids, and consider 
the Q-tangent prolongation of i.e b: T e E — > T E E 1 . Let and £' v , S 
and S', and A and A' , be the vertical liftings, the vertical endomorphisms, and the 
Liouville sections on E and E', respectively. Then, 

(1) T*$(^(a, b )) = £' v ($(a), $(&)), for all ( a,b)£Ex M E, 

(2) T $ <f> oA = A'o$, 

(3) T 4 $oS = S' oT $ $. 


Proof. For the first property, we notice that both terms are vertical, so that we 
just have to show that their action on functions coincide. For every function f £ 
C°°(E'), we deduce that 


p ,1 (T*$(f v (a, &)))/' 


T3>{p\e (a, b)))f = T$(b v a )f' = b v a {f a $) 


-f'($(a + tb)) =-/ , (4>(a) + t4>(6)) 
at t =o at 


t=0 


For the second property, we have A (a) = (a, a) so that applying the first property 

it follows that 


T*$(A(a)) = T*<h(r (a, a)) = r (<&(«), <&(a)) = A'(<b(a)). 
Finally, for any z = (a, b, V) £ T e E , we obtain that 

T*$(S(*)) = T*$(r (a, 6)) = r ($(«),$(&)) 

= S' (4>(a), $(6), T$(P)) = S' (T $ $(z)), 
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which concludes the proof. □ 

Proposition 4.2. Let L £ C°°(E) be a Lagrangian function, 9 l the Cartan form 
and u>l = —dOr- Let •!>: E —> E' be a Lie algebroid morphism and suppose that 
L = L' o <f>, with L' £ C°°(E ') a Lagrangian function. Then, we have 

(1) ( T*®ye L , = e L , 

(2) ( t *$)* w l , = u L , 

(3) (T*<f>)*E L , = e l , 

( 4 ) G $(o)( $ ( 6 )> $ ( c )) = G a ( b > c), /or every and every b,c& E T ^ a) . 

Proof. Indeed, for every Z £ T E E we have 

((T*$)*0 L ,,Z) = (0 L ,,T*$(Z)) = (dZ/,S'(T*$(Z))> = (dL',T*$(S(Z))) 

= ((T*$)*dL', S(Z)) = (d(T*$)*L', S{Z )) = (d(L' o $), S(Z)) 
= (dL,s(z)) = (e L ,z), 

where we have used the transformation rule for the vertical endomorphism. The 
second property follows from the fact that T 3 ’<I> is a morphism, so that (T $ <E>)*d = 
d(T $ <f>)*. The third one follows similarly and the fourth is a consequence of the 
second property and the definitions of the tensors G L and G L '. □ 

Let T be a SODE and L £ C°°(E ) be a Lagrangian. For convenience, we define 
the 1-form SrL £ Sec ((T e E)*) by 

(S r L, Z) = (dE L - i T u L , Z) = ( dE L , Z) - u L (T, Z), 

for every section Z of T e E. We notice that T is the solution of the free dynamics 
if and only if SrL = 0. On the other hand, notice that the 1-form SrL is semibasic, 
because T is a sode. 

Proposition 4.3. Let T be a SODE in E and T' a SODE in E'. Let L £ C°°(E) 
and L' £ C 00 (E') be Lagrangian functions defined on E and E', respectively, such 
that L = L' o$. Then, 

(6 r L~- (T*$)*6r>L',Z) =w L -(r'o$ —T 4 $oT,T $ $(Z)), (4.1) 

for every section Z ofT E E. 

Proof. Indeed, from (T^i>)*dE^ = dEf, we have that 

(S r L - (T*$)*Sr>L', Z) = ((T*$)*i r /Wj/ - »twl, Z) 

= <(T*$)*ir'Wi/ -i r (T*$)*Wi/,Z) 

= wi,/(r' o $ - o T, r*$(Z)), 

which concludes the proof. □ 

4.1. Reduction of the free dynamics. Here, we build on Propositions l4.2l and l4.3l 
to identify conditions under which the dynamics can be reduced under a morphism 
of Lie algebroids. We first notice that, from Proposition l4.21 if $ is fiberwise surjec¬ 
tive morphism and L is a regular Lagrangian on E, then L' is a regular Lagrangian 
on E' (note that : T e E -> T E E' is a fiberwise surjective morphism). Thus, 
the dynamics of both systems is uniquely defined. 

Theorem 4.4 (Reduction of the free dynamics). Suppose that the Lagrangian func¬ 
tions L and L' are 4>- related, that is, L = L' o4>. If <f> is a fiberwise surjective mor¬ 
phism and L is a regular Lagrangian then L' is also a regular Lagrangian. Moreover, 
ifTr and Tl' are the solutions of the free dynamics defined by L and L' then 

T^of^f^of. 
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Therefore, if aft) is a solution of the free dynamics defined by L, then 4>(a(f)) is a 
solution of the free dynamics defined by L'. 

Proof. If Fl and T£/ are the solutions of the dynamics, then dr L L = 0 and Sr L , L' = 

0 so that the left-hand side in equation SH) vanishes. Thus 

u L ,(r L , o$-T $ $or Ll r $ $(z)) = o, 

for every Z £ Sec (T e E). Therefore, using that L' is regular and the fact that 
is a fiberwise surjective morphism, we conclude the result. □ 

We will say that the unconstrained dynamics T rj is the reduction of the un¬ 
constrained dynamics I’/, by the morphism $. 

4.2. Reduction of the constrained dynamics. The above results about reduc¬ 
tion of unconstrained Lagrangian systems can be easily generalized to nonholonomic 
constrained Lagrangian systems whenever the constraints of one system are mapped 
by the morphism to the constraints of the second system. Let us elaborate on this. 

Let (L,D) be a constrained Lagrangian system on a Lie algebroid E and let 
(. L',D') be another constrained Lagrangian system on a second Lie algebroid E'. 
Along this section, we assume that there is a fiberwise surjective morphism of Lie 
algebroids 4>: E —> E' such that L = L' o $ and <1>(D) = D'. The latter condition 
implies that the base map is also surjective, so that we will assume that $ is an 
epimorphism (i.e., in addition to being fiberwise surjective, the base map p is a 
submersion). 

As a first consequence, we have ($(&), $(c)) = G^’ D (b, c), for every a £ D 

and every b,c £ D 7r ( Q ), and therefore, if (L, D) is regular, then so is (L', D'). 

Lemma 4.5. With respect to the decompositions T e E\b = T e D(BF and T E E'\e>i = 
T E ' D' ®F', we have the following properties: 

(1) T*${T e D) = T e 'D’, 

(2) T^^(F) = F', 

(3) If P,Q and P',Q' are the projectors associated with ( L,D ) and ( L',D '), 

respectively, then P' o = T $ <f> o P and Q' o o Q. 

With respect to the decompositions T e E\d = T D D © (T D D) E and T E E’\yo' = 
T D D' © (T D D') e we have the following properties: 

(4) T®<&{T d D) = T d 'D’, 

(5) T' s> ®((T d D) ± ) = (T D 'D') e , 

(6) If P,Q and P',Q' are the projectors associated with ( L,D) and ( L',D'), 

respectively, then P' o o P and Q' o o Q. 

Proof. From the definition of T $ 4>, it follows that 

(T' l ‘<I>)(T E D) C T E 'D', (T $ $)(T I) £i) C T D 'D'. 

Thus, one may consider the vector bundle morphisms 

r $ 4> : T e D T E 'D\ T $ 4> : T°D -> T°'D'. 

Moreover, using that is fiberwise surjective and that <p is a submersion, we deduce 
that the rank of the above morphisms is maximum. This proves (1) and (4). 

The proof of (5) is as follows. For every a' £ D ', one can choose a £ D such that 
<I>(a) = a', and one can write any element w' £ T E D' as w' = T $< F(m;) for some 
w £ T E D. Thus, if z £ (T E D) ± , for every w' £ T E 'D' we have 

co L '{r^^z),w') = =w L {z,w) = 0, 
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from where it follows that £ (T D D')^. In a similar way, using that 

T $ 4> : (T e E) |£) —► ( T e E') j£>/ is fiberwise surjective, (2) in Proposition 14.21 and 
(4), we obtain that {T D 'D') E Q {T‘ s, ^){{T d D) ± ). 

For the proof of (2) we have that 

T $ 4>(F) = n Vei(T E E)) C (T D 'D') E n \ei(T E 'E') = F'. 

Thus, using that T $ 4> : (T E F)|£) —> {T E 'E')\ D > is fiberwise surjective, the fact 
that (T e E) | d = T e D ® F and (1), it follows that 

(T E 'E')\ D , = T E> D' © (T*$)(F). 

Therefore, since (T B E')^i = T E D' © F', we conclude that (2) holds. 

Finally, (3) is an immediate consequence of (1) and (2), and similarly, (6) is an 
immediate consequence of (4) and (5). □ 

From the properties above, we get the following result. 

Theorem 4.6 (Reduction of the constrained dynamics). Let (L,D) be a regular 
constrained Lagrangian system on a Lie algebroid E and let (L ', D') be a constrained 
Lagrangian system on a second Lie algebroid E'. Assume that a fiberwise surjective 
morphism of Lie algebroids 4>: F —> F' exists such that L = L' o$ and 4 ’(D) = 
D'. If is the constrained dynamics for L and is the constrained 

dynamics for L', then T*4> oV^ LD \ = r°4>. If a(t) is a solution of Lagrange- 
d’Alembert differential equations for L, then 4>(a(f)) is a solution of Lagrange- 
d’Alembert differential equations for L'. 

Proof. For the free dynamics, we have that T $ 4> oTi = T l' o <f>. Moreover, from 
property (3) in Lemma 14.51 for every a £ D, we have that 

T*$(r {LtD) (a)) = T*<t>(P(T L (a))) = P'(T**(T L (a))) 

= P'(r i ,($(a)))=r (i ,, D0 ($(o)), 

which concludes the proof. □ 

We will say that the constrained dynamics F (l'.d 1 ) is the reduction of the 
constrained dynamics T( L D ) by the morphism 4>. 

Theorem 4.7. Under the same hypotheses as in Theorem El we have that 
{f'o^g l o^} nh = {f',g'y nh o^, 

for f',g' £ C°°(D'), where {-,-} n h (respectively, is the nonholonomic 

bracket for the constrained system ( L,D ) (respectively, ( L’,D’)). In other words, 
4> : D —> D' is an almost-Poisson morphism. 

Proof. Using (2) in Proposition l4.2l and the fact that $ is a Lie algebroid morphism, 
we deduce that 

(ix flo t(T*<f>)*u L >) = i Xf ,u>L> o 4-. 

Thus, since T $ 4> is fiberwise surjective, we obtain that 

T* 4> o X/- 04 . = X r o <P. 

Now, from (13.111) and Lemma B~5l we conclude that 

{f'o<S>,g'o$} nh = {f',g'}' nh o<S>. 

□ 

One of the most important cases in the theory of reduction is the case of reduction 
by a symmetry group. In this respect, we have the following result. 
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Theorem 4.8 ([131131]). Let q^,: Q —> M be a principal G-bundle, letr: E —> Q be 
a Lie algebroid, and assume that we have an action of G on E such that the quotient 
vector bundle E/G is well-defined. If the set Sec (E) G of equivariant sections of E 
is a Lie subalgebra ofSec(E), then the quotient E' = E/G has a canonical Lie 
algebroid structure over M such that the canonical projection q^: E —v E/G, given 
by a i—> [o\g, is a (fiberwise bijective) Lie algebroid morphism over q q. 

As a concrete example of application of the above theorem, we have the well- 
known case of the Atiyah or Gauge algebroid. In this case, the Lie algebroid E is 
the standard Lie algebroid TQ —> Q, the action is by tangent maps gv = Ti[> g (v), 
the reduction is the Atiyah Lie algebroid TQ/G —> Q/G and the quotient map 
q^: TQ —> TQ/G is a Lie algebroid epimorphism. It follows that if L is a G- 
invariant regular Lagrangian on TQ then the unconstrained dynamics for L projects 
to the unconstrained dynamics for the reduced Lagrangian L'. Moreover, if the 
constraints D are also G-invariant, then the constrained dynamics for (L, D) reduces 
to the constrained dynamics for ( L',D/G ). 

On a final note, we mention that the pullback of the distributional equation 
i T 'OJ L '' D ' - e L '’ D ' = 0 by 7^$ is precisely (i T u L ' D - e L ' D ) o = 0. 


4.3. Reduction by stages. As a direct consequence of the results exposed above, 
one can obtain a theory of reduction by stages. In Poisson geometry, reduction by 
stages is a straightforward procedure. Given the fact that the Lagrangian counter¬ 
part of Poisson reduction is Lagrangian reduction, it is not strange that reduction 
by stages in our framework becomes also straightforward. 

The Lagrangian theory of reduction by stages is a consequence of the following 
basic observation: 


Let 4>i: Eq —> E\ and <f >2 '■ E\ —> E 2 be a fiberwise surjective mor- 
phisms of Lie algebroids and let <f>: Eo —:► E 2 be the composition 
4> = $2 o $ 1 . The reduction of a Lagrangian system in Eo by $ can 
be obtained by first reducing by <f>i and then reducing the resulting 
Lagrangian system by <f >2 • 

This result follows using that T $ 4> = T $2 <f >2 o T $1< f> 1 . Based on this fact, one 
can analyze one of the most interesting cases of reduction: the reduction by the 
action of a symmetry group. We consider a group G acting on a manifold Q and a 
closed normal subgroup N of G. The process of reduction by stages is illustrated 
in the following diagram 


Total reduction 



T Q '■ Eq = TQ —► M 0 = Q 


First reduction 


■/N 


t 1 :E 1 = TQ/N -► M 1 = Q/N 


Second reduction 


•/( G/N ) 


r 2 : E 2 = (TQ/N)/(G/N) - M 2 = (Q/N)/(G/N) 


In order to prove our results about reduction by stages, we have to prove that 
Eq, Ei and E 2 are Lie algebroids, that the quotient maps <I>i: Eq —■> E\, $2 : E\ —> 
E 2 and 4> : Eq — > E 2 are Lie algebroids morphisms and that the composition 4>i o4> 2 
equals to 4>. Our proof is based on the following well-known result (see ns). which 
contains most of the ingredients in the theory of reduction by stages. 

Theorem 4.9. 7[15| ) Let q®: Q —> M be a principal G-bundle and N a closed 
normal subgroup of G. Then, 

(1) q^: Q —> Q/N is a principal N-bundle, 
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(2) G/N acts on Q/N by the rule [g]jv[g]jv = [9<i]n, 

(3) cIq/%'- Q/N —> (Q/N) / (G/N) is a principal (G/N)-bundle. 

(4) The map i: Q/G —> (Q/N)/(G/N) defined by [g]c i—> [[^jvJg/jv a diffeo- 

morphism. 

Building on the previous results, one can deduce the following theorem, which 
states that the reduction of a Lie algebroid can be done by stages. 

Theorem 4.10. Let : Q —> M be a principal G-bundle and N be a closed normal 
subgroup of G. Then, 

(1) t T q /g : TQ/G —> Q/G is a Lie algebroid and qJP: TQ —> TQ/G is a Lie 
algebroid epimorphism, 

(2) ttq/n'- TQ/N —> Q/N is a Lie algebroid and qj^: TQ —> TQ/N is a Lie 
algebroid epimorphism, 

(3) G/N acts on TQ/N by the rule [<7]jv[V]iv = \gv\n, 

(4) T(tq/ N )/(g/ N ) '■ (TQ/N)/(G/N) —> (Q/N)/(G/N) is a Lie algebroid and 

'• TQ/N —i► (TQ/N)/(G/N) is a Lie algebroid epimorphism, 

(5) The map I: TQ/G — ► (TQ/N)/(G/N) defined by [u]g i—► [[u]jv]g/jv an 
isomorphism of Lie algebroids over the map i. 

Proof. The vector bundle t T q/q ■ TQ/G —► Q/G (respectively, ttq/n '■ TQ/N — ► 
Q/N) is the Atiyah algebroid for the principal G-bundle q^ : Q —» Q/G (respec¬ 
tively, q® : Q —> Q/N), so that (1) and (2) are obvious. Condition (3) is just 
condition (2) of Theorem 14.91 applied to the principal TV-bundle TQ —> TQ/N. To 
prove condition (4), we notice that the action of G/N on the Lie algebroid TQ/N 
is free and satisfies the conditions of Theorem 14.81 Finally, the Lie algebroid mor¬ 
phism j: TQ —► TQ/N is equivariant with respect to the G-action on TQ and 
the (G/^V)-action on TQ/N. Thus, it induces a morphism of Lie algebroids in the 
quotient. It is an isomorphism since it is a diffeomorphism by Theorem 14.91 □ 


The following diagram illustrates the above situation: 



$ 2 


(TQ/N)/(G/N) 


a Q/N 

Vg/N 


(Q/N)/ (G/N) 


In particular, for the unconstrained case one has the following result. 

Theorem 4.11 (Reduction by stages of the free dynamics). Let q®: Q —> Q/G be 
a principal G-bundle, and N a closed normal subgroup of G. Let L be a Lagrangian 
function on Q which is G-invariant. Then the reduction by the symmetry group G 
can be performed in two stages: 

1. reduce by the normal subgroup N, 

2. reduce the resulting dynamics from 1. by the residual symmetry group G/N. 


Since the dynamics of a constrained system is obtained by projection of the free 
dynamics, we also the following result. 

Theorem 4.12 (Reduction by stages of the constrained dynamics). Let q®: Q —> 
Q/G be a principal G-bundle and N a closed normal subgroup of G. Let (L,D) be 
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a G-invariant constrained Lagrangian system. Then the reduction by the symmetry 
group G can be performed in two stages: 

1. reduce by the normal subgroup N, 

2. reduce the resulting dynamics from 1. by the residual symmetry group G/N. 


5. The momentum equation 

In this section, we introduce the momentum map for a constrained system on a 
Lie algebroid, and examine its evolution along the dynamics. This gives rise to the 
so-called momentum equation. 


5.1. Unconstrained case. Let us start by discussing the unconstrained case. Let 
te : E —► M be a Lie algebroid over a manifold M and L : E —> R be a regular 
Lagrangian function. Suppose that tk '■ K —» M is a vector bundle over M and that 
T : K —> E is a vector bundle morphism (over the identity of M) between K and 
E. Then, we can define the unconstrained momentum map J(ln>) : E —> K* 
associated with L and T as follows 


J(L,it)(a) e K*, for a G E x , 

and 

= 4 L(a + m(k)) = for k G K x . 

dt\t=o 

If a : M —> K is a section of tk '■ K —> M then, using the momentum map 
we may introduce the real function J* L ^ : E — >■ R given by 

= J lL^){a)(a{x)) = y(cr(x))l{L), for a G E x . (5.1) 

Theorem 5.1 (The unconstrained momentum equation). Let T l be the Euler- 
Lagrange section associated with the regular Lagrangian function L : E —> R. If 
a : M —> K is a section of tk '■ K —> M and (T o a) c G Sec(T E E) is the complete 
lift of (tE 1 o a) G Sec(E), we have that 

< d TEE J^ ) ,T L >=< rt,($o<7) c >, (5.2) 

where d T E is the differential of Lie algebroid T e E —> E. In particular, if 
< d E E L, (T o cr) c >= 0, then the real function J° L ^ is a constant of the mo¬ 
tion for the Lagrangian dynamics associated with the Lagrangian function L. 

Proof. Let S : T e E —> T e E be the vertical endomorphism. If (Tocr)" G Sec(T E E) 
is the vertical lift of (fo<r) G Sec(E) then, using (15.ID and the fact that 5(H'ocr) c = 
(T o a) v , it follows that 

j(L,*)=0L«*°*n (5-3) 

where 0^ is the Cartan 1-form associated with L. 

Thus, from (|5.3D . we deduce that 

d TEE J[ L ,*) = E „y e L + *(H/ o a) c (w L ), 

ujl being the Cartan 2-form associated with L. 

Therefore, if Ee ■ E —> M is the Lagrangian energy, we obtain that 

<d rBE j( w T L > = <d TEE (e L (r L )),^oay>^<d TEE E L ,(*oay> 

-<e L ,[^o a y,T L ]>. 


Now, from (12.111) and since Tl is a SODE section, it follows that 

©L(r L ) =< d rEE L, a >, < e L , [(T o cr) c , r L ] >= o, 


(5.4) 
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where A G Sec{T E E) is the Liouville section. Consequently, using (IQ1) we deduce 
that (15.21) holds. □ 

Remark 5.2 (Conservation of momentum on TM). Let L : TM —> R be an 
standard regular Lagrangian function on T M. Suppose that G is a Lie group with 
Lie algebra g and that ^:GxM-»Misa (left) action of G on M. Then, we may 
consider the trivial vector bundle over M 

K = M x g —> M 

and the vector bundle morphism T : K —► TM (over the identity of M) defined by 

= £m(z), (5.5) 

where G X(M) is the infinitesimal generator of the action ip associated with 

£ G 0- 

A direct computation proves that the (unconstrained) momentum map Jl,^) '■ 
E = TM —> K* = Mxg* associated with L and is given by 

J(l, 9 )(vx) = (x,J(v x )), for v x G T X M , 

where J : TM —> g* is the standard momentum map associated with L and the 
action ip defined by 

J{ v x){& = 4 L ( v * + t£ M {x)), for v x G T X M and £ G g 
at |t=o 

(see, for instance, DP)- 

Now, each £ G g dehnes a (constant) section a of the vector bundle K = M x g —> 
M and the real function J? L is just the momentum : TM —> R in the direction 
of C 

On the other hand, if r/ G g, then the infinitesimal generator T]tm of the tangent 
action Tip : G xTM ^ TM associated with 77 is the (standard) complete lift 77 ^ G 
X(TM) of t]m- Therefore, using Theorem 15.11 we deduce a well-known result [TJ: 
“If the Lagrangian function L : TM —> K. is invariant under the tangent action Tip 
of G on TM then, for every ( G g, the momentum : TM —> R in the direction 
of £ is a constant of the motion of the Lagrangian dynamics.” • 

5.2. Constrained case. Next, let us discuss the constrained case. Suppose that 
L : E —> R is a regular Lagrangian function on a Lie algebroid te '■ E —> M, that 
tk : K M is a vector bundle over M and that T : K —> E is a vector bundle 
morphism (over the identity of M) between K and E. 

In addition, let td '■ D —> M be a vector subbundle of te ■ E —> M such that 
the nonholonomic Lagrangian system (L, D) is regular. 

If x is point of M we consider the vector subspace K E of K x given by 

K x = {k G K x /*{k) G D x }. 

We will denote by i x : K x —7 K x the canonical inclusion, by i* : K* —> ( K x )* the 
canonical projection and by K D and ( K D )* the sets 

k d = u (K D y = u (K! r y. 

xGM xEM 

Then, we define the nonholonomic momentum map J(l,d , t ) : E —> (K D )* 
associated with the system (L , D ) and the morphism if as follows 

{J(L,D,*))\E„ =i*x° for X G M. 

Now, if a : M — > K is a section of tk '■ K —> M such that cr(x) G K x , for all 
x G M, we may introduce the real function J° L : E —> R given by 

J (L,D,*)( a ) = J (L,D,^){a)(cr(x)) : for a G E x , 
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that is, J( LD 

Theorem 5.3 (The nonholonomic momentum equation). LetT m be the solution 
of the constrained dynamics for the nonholonomic Lagrangian system ( L,D ). If 
a : M —> K is a section of tk : K —> M such that &{x) £ K E , for all x £ M, and 
('I' o a) c £ Sec(T E E) is the complete lift of ('I 1 o a) £ Sec(E) then we have that 

< d T ^(l,d) >=< d T e L, (\& o a) c >|£,, (5.6) 

where d T D (respectively, d T E ) is the differential of Lie algebroid T E D —> D 
(respectively, T E E —> E). In particular, if < d T E L, (if o cr) c >\d = 0, then the 
real function J( L D ^ is a constant of the motion for the constrained dynamics 
associated with the nonholonomic Lagrangian system ( L,D ). 

Proof. Denote by j : D —> E and by J : T E D —> T E E the canonical inclusions and 
by Q : T E E —> F the corresponding projector, where F = w£ 1 (I) 0 ) (see Section 
USD - Then, as we know, 

r (l,d) = (Tl — QLl)\d- 

Moreover, the pair {J,j) is a Lie algebroid monomorphism which implies that 

d T D = {J,j)*(d T E 

Thus, using that J( L D ^ = J( L ^ and proceedings as in the proof of Theorem 15.II 
we deduce that 

<d TED ((J( w) ) lD ), r (i>D) > = <d rEE L,(^oaY> ]D 

+ (i(* O anco L ))(QT L )} lD . 

Now, since S(QTl) = 0, then 5[(^ o <j) c ,QTl] = 0 (see <|2.12|> ) and it follows that 
9l(QT l ) = 0, e L [^oa) c ,QT L ] = 0. 

Therefore, 

(C^ E r 9 L )(QT L ) = 0. (5.8) 

On the other hand, we have that 

(i{QT l)u l )\d = S*(at( L ' D )), with £ Sec((T E D)°). 

Consequently, 

{(i{ip o a) c u L )(QT L )} lD = o cr)J D ). 

But, since \F o a is a section of : D —> M, it follows that (>F o a)( D is a section 
of T E D —> D. This implies that 

{(i{^ oa) c Lo L )(QT L )} lD = 0. (5.9) 

Finally, using (15.71) . (15.81) and (15.91) . we conclude that (15.61) holds. □ 

Remark 5.4 (Nonholonomic momentum equation on TM and horizontal symme¬ 
tries). Suppose that L : TM —> R is an standard regular Lagrangian function on 
E = TM and that : G x M —> M is a (left) action of a Lie group G on M. Then, 
we consider the trivial vector bundle tk '■ K = M x g —> M and the vector bundle 
morphism *F : K —> TM (over the identity of M) defined by (15.51) . 

Now, let I? be a vector subbundle (over M) of the vector bundle tm '■ TM —> M, 
that is, D is a distribution on M, and assume that the nonholonomic Lagrangian 
system (L, D) is regular. If a; is a point of M, we have that K E = {&} x g x , where 
0 * is the vector subspace of g given by 

Q x = {£ e b/€m{x) £ D x }. 
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We also remark that the sets K D and ( K D )* may be identified with the sets 

g D = 1J 0 "’ (s D )* = U (0 X )*- 

iEM i6« 

Under this identification, the nonholonomic momentum map J(l.d,< b) : E —► ( K D )* 
associated with the system ( L,D ) and the morphism T is just the standard non¬ 
holonomic momentum map J nh : TM —► (g D )* associated with the system (L,D) 
and the action tjj (see 000). 

Now, if £ : M —► g is an smooth map the £ defines, in a natural way, a section 

C7 


• 

J~ h 


M ^ K = M x g of the vector bundle tk '■ K = M x g —> M. We denote by 

CTr- 

the real function Dt j : E —> R and by Sj- the vector field U/oct^- 


on M. Then, using Theorem 15.31 we deduce a well-known result (see 000): “If 
F (l,d) is the solution of the constrained dynamics for the nonholonomic system 
(L 1 D ), we have that 

TL,D((Jf )\d) = (3 f )fb(L).” 

The above equality is an intrinsic expression of the standard nonholonomic 
momentum equation. In addition, using again Theorem 15.31 we also deduce 
another well-known result (see 000): “If the Lagrangian function I : TM-> R 
is invariant under the tangent action of G on TM and ( 6 g is a horizontal 
symmetry (that is, £ £ g x , for all x £ M ) then the real function (J~ h )\D is a 

constant of the motion for the constrained Lagrangian dynamics, where £ : M —> g 
is the constant map 

£(x) = £, for all x G M.” • 

6. Examples 

As in the unconstrained case, constrained Lagrangian systems on Lie algebroids 
appear frequently. We show some examples next. 

6.1. Nonholonomic Lagrangian systems on Lie algebras. Let g be a real 
algebra of finite dimension. Then, it is clear that g is a Lie algebroid over a single 
point. Now, suppose that (Z, D) is a nonholonomic Lagrangian system on g, that is, 
l : g —> R is a Lagrangian function and 0 is a vector subspace of g. If w : I —> g is 
a curve on g then 

dl(w(t)) G = g*, Vt G I, 

and thus, the map dl o to may be considered as a curve on g* 

dl o uj : I —> g*. 

Therefore, 

(dl o uj)'{t) G T dl{uW) g* ^ g*, Vt G I. 

Moreover, from EZD, it follows that to is a solution of the Lagrange-d’Alembert 
equations for the system (Z, t>) if and only if 

(dl o tL>)'(t) — ad^,^(dl(uj(t))) G D°, u(t) G 0, Vt (6-1) 

where ad* : g x g* —» g* is the infinitesimal coadjoint action. 

The above equations are just the so-called Euler-Poincare -Suslov equations 
for the system (t,0) (see [2SJ). We remark that in the particular case when the 
system is unconstrained, that is, 0 = g, then one recovers the the standard Euler- 
Poincare equations for the Lagrangian function l : g —> R. 

If G is a Lie group with Lie algebra g then nonholonomic Lagrangian systems on 
g may be obtained (by reduction) from nonholonomic LL mechanical systems with 
configuration space the Lie group G. 
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In fact, let e be the identity element of G and I : g —> g* be a symmetric positive 
definite inertia operator. Denote by g e : g x g —> R the corresponding scalar product 
on g given by 

g e (ui,ui') =< I(w),w' >, for w,u'eg = T e G. 

g e induces a left-invariant Riemannian metric g on G. Thus, we way consider the 
Lagrangian function L : TG —> R defined by 

L{v h ) = ^gh(v h ,v h ), for Vh S T h G. 

In other words, L is the kinetic energy associated with the Riemannian metric g. 

Now, let D be a left-invariant distribution on G. Then, since L is a left-invariant 
function, the pair (L, D) is an standard nonholonomic LL system in the terminology 

of [25] , 

On the other hand, the Lagrangian momentum map $ : TG —> g given by 
${vh) = {T h l h -t){yh), for v h G T h G 

is a fiberwise bijective morphism of Lie algebroids. Moreover, if l = L g and 0 = D e 
then the pair (Z,5) is a nonholonomic Lagrangian system on g and 

l o <f> = L and $(£>) = 5. 

Thus, the system (Z,5) is regular. In addition, if v : I —► TG is a solution of 
the Lagrange-d’Alembert equations for the system (L, D) then, using Theorem 14.61 
we deduce that the curve <I> o v : I —> g is a solution of the Lagrange-d’Alembert 
equations for the system (Z,S). 

We remark that 

l (u) = ^. 9 e(w,w) = i < I(w),w >, for w G g. 

Therefore, if ui : I —> g is a curve on g, we have that 

[dl o uj)(t) = I(w(t)), for all t 

and, using (EH), it follows that w is a solution of the Lagrange-d’Alembert equations 
for the system (7,5) if and only if 

uj — I^ 1 (ad* ( t )I(w(f))) G 5^, uj{t) G 5, for all t, 

where 5^ is the orthogonal complement of the subspace 0, that is, 

5 X = {li/ Gg/< I(w'),w >= 0,Vw G 5}. 

Two simple examples of the above general situation are the following ones. 

The Suslov system. The most natural example of LL system is the nonholo¬ 
nomic Suslov problem , which describes the motion of a rigid body about a fixed 
point under the action of the following nonholonomic constraint: the body angular 
velocity vector is orthogonal to a some fixed direction in the body frame. 

The configuration space of the problem is the group G = 50(3). Thus, in this 
case, the Lie algebra g may be identified with M 3 and, under this identification, the 
Lie bracket on g is just the cross product x on R 3 . 

Moreover, if I : R 3 —> (R 3 )* = R 3 is the inertia tensor of the body then a curve 
u> : I —> R 3 on R 3 is a solution of the Euler-Poincare-Suslov equations for the 
system if and only if 

w = I" 1 ((Iw) X w)) + AI- 1 (T), <w,T>=0, (6.2) 

where A is the Lagrange multiplier, T is a fixed unit vector in R 3 and < *, * > is 
the standard scalar product in R 3 . Since the nonholonomic system is regular, the 
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Lagrange multiplier A is uniquely determined. In fact, differentiating the equation 
< u, T >= 0, we find 

<iuxu, i _1 r > 

< r,i-!r > 

and, consequently, Eqs. (16.21) are equivalent to 

uj = i _1 (< iw,r > u> x i - 1 r), < w,r >= o. 

Multidimensional generalizations of the Suslov problem have been discussed by 
several authors (see EH HflES]). 


The Chaplygin sleigh. The Chaplygin sleigh is a rigid body sliding on a hor¬ 
izontal plane. The body is supported at three points, two of which slide freely 
without friction while the third is a knife edge, a constraint that allows no motion 
orthogonal to this edge. This mechanical system was introduced and studied in 
1911 by Chaplygin [TB] (see also [57]). 

The configuration space of this system is the group SE( 2) of Euclidean motions 
of the two-dimensional plane R 2 . As we know, we may choose local coordinates 
(0, x, y) on SE( 2). 9 and (x, y) are the angular orientation of the blade and position 
of the contact point of the blade on the plane, respectively. 

Now, we introduce a coordinate system called the body frame by placing the 
origin at the contact point and choosing the first coordinate axis in the direction 
of the knife edge. Denote the angular velocity of the body by w = 9 and the 
components of the linear velocity of the contact point relative to the body frame by 
V \, v 2 ■ The set (u>,vi,v 2 ) is regarded as an element of the Lie algebra se(2). Note 
that 


V\ = x cos 9 + y sin (9, V 2 = y cos 9 — x sin 9. 


The position of the center of mass is specified by the coordinates (a, b) relative to 
the body frame. Let m and J denote the mass and moment of inertia of the sleigh 
relative to the contact point. Then, the corresponding symmetric positive definite 
inertia operator I : se(2) —> se(2)* and the reduced nonholonomic Lagrangian 
system (Z,0) on se(2) are given by 

( J + m(a 2 + b 2 ) —bm am \ / u> \ 

—bm m 0 j ( z;i ] , 

am 0 m J \ V 2 J 

\[(J + m(a 2 + b 2 ))u} 2 + m(v\ + v 2 ) 

—2mbujv\ + 2 amujV 2 \ 1 

5 = {(w,i>i,u 2 ) € se(2)/u 2 = 0}, 


l(u>,v 1 ,v 2 ) = 
l(u,vi,v 2 ) = 


(see [IE]). Thus, the Lagrange-d’Alembert equations for the system (Z,0) are 


<jj 


Vi 

V 2 


amuj 

J + ma 2 
au> 

J + ma 2 
0. 


( bco - ui), 

((J + to(o 2 + b 2 ))uj 


mbv i), 


Multidimensional generalizations of the Chaplygin sleigh were discussed in [2B] (see 
also m and ESI). 





NONHOLONOMIC LAGRANGIAN SYSTEMS ON LIE ALGEBROIDS 


31 


6.2. Nonholonomic LR systems and right action Lie algebroids. Here, we 
show how the reduction of a nonholonomic LR system produces a nonholonomic 
Lagrangian system on a right action Lie algebroid. 

Let us start by recalling the definition of a right action Lie algebroid (see [32]). 
Let (F, [•, ■] jr, pf) be a Lie algebroid over a manifold N and 7r : M —> N be a smooth 
map. A right action of F on n : M —> TV is a R-linear map 

: Sec(F) -► X{M), X G Sec(F) -> V(X) G X(M) 

such that 

*(fX) = (/ ° tt)®(X), ([X, Y\ f ) = [*(*), (Y)], 

(T m 7r)(^(X)(m)) = p F (X(Tr(m))), 

for / G C°°(N), X, Y G Sec(F) and m G M. If T : Sec(E) —► X(M) is a right 
action of F on 7r : M —> TV and te : F —> TV is the vector bundle projection then 
the pullback vector bundle of F over 7r, 

E = n*F = {(to, /) G M X F/r F (f) = 7r(m)} 

is a Lie algebroid over M with Lie algebroid structure ([• i -]e , Pe) which is charac¬ 
terized by 

[X,Y] e = [X,Y] f ott, p E (X)(m) = tf(X)(m), 
for I,Y G Sec(E) and to G M. The triple (F, [•, -]s, Pe) is called the right action 
Lie algebroid of F over n and it is denoted by 7 t^F (see [32]). 

Note that if the Lie algebroid F is a real Lie algebra g of finite dimension and 
7r : XL —> {a point} is the constant map then a right action of g on 7r is just 
a right infinitesimal action : g —> X(M) ofg on the manifold M. In such 
a case, the corresponding right action Lie algebroid is the trivial vector bundle 
pr\ :Mxg^ M. 

Next we recall the definition of a nonholonomic LR system following [2E, .35]. 
Let G be a compact connected Lie group with Lie algebra g and < •, • >: g x g —> R 
be an Hd^-invariant scalar product on g. Now, suppose that X : g —*• g is a inertia 
operator which is symmetric and definite positive with respect to the scalar product 
< ■, ■ >. Denote by g the left-invariant Riemannian metric given by 

9ii{v h ,v h >) =<l(T h l h -i(v h )),(T h l h -i)(v' h ) > (6.3) 

for h G G and Vh, v' h G T^G. 

Then, the Lagrangian function L : TG —> R of the system is 

L(vh) = t -; 9 h(vh,v h ) - V(h), for v h G T h G , (6.4) 

V : G —> R being the potential energy. The constraint distribution D is a right- 
invariant distribution on G. Thus, if e is the identity element of G and c) = D ei we 
have that 

D h = ( T e r h )(* ) = (T e l h )(Ad h - i(D)), for h G G (6.5) 

where Ad : G x g —> g is the adjoint action. 

The nonholonomic Lagrangian system ( L,D ) on TG is called a nonholonomic 
LR system in the terminology of [20 [35]. Note that, since L is a Lagrangian 
function of mechanical type, the system (L, D) is regular. Now, assume that 

5 = 0^= {a/ G g/ < ui, to 1 >= 0, Vw G 0} 

is a Lie subalgebra of g, that S' is a closed Lie subgroup of G with Lie algebra s 
and that the potential energy V is S-invariant. 

Next, let us show that the nonholonomic LR system (L, D) may be reduced to a 
nonholonomic Lagrangian system on a right action Lie algebroid. In fact, consider 
the Riemannian homogeneous space M = S \ G and the standard transitive right 
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action ip of G on M = S\G. Denote by 'J/ : g —> X(S \ G) the corresponding right 
infinitesimal action of g on S \ G. Then, tl/ induces a Lie algebroid structure on the 
trivial vector bundle pri :S\Gxg^>S\G. 

On the other hand, using that the potential energy V is S invariant, we deduce 
that V induces a real function V : S\G —on S\G such that 

b o 7 T = V, (6.6) 

where 7r : G —> S \ G is the canonical projection. Thus, we can introduce the 
Lagrangian function L:S\G x g —> Ron the action Lie algebroid pr± : S\G x g —> 
S\G defined by 

L{h,u) = i < I(lo), u> > — V(h), for h £ S \ G and oj £ g. (6.7) 

Now, for every h £ G, we consider the subspace d(h) of g given by 

d(h) = Adh~i (0). (6-8) 

The dimension of d(h) is equal to the dimension of 0 . Moreover, since < •, • > is 

Ad^-invariant, if follows that 

0 (/i) = {Ad h -i (s)) 1 - = {a / eg/ < a/, Ad h -1 (w) >= 0, Vw £ s}. 

In particular, we have that 

h(s) = s 1 ' = h, Vs S S 

which implies that 5(s/i) = D(h), for all h £ G. 

Therefore, we can define a vector subbundle D of the Lie algebroid pr\ : S\G x 
0 —> S \ G as follows 

= {h} x V(h), for h G S\G (6.9) 

with h £ G and 7r (h) = h. Consequently, the pair (L, D) is a nonholonomic La¬ 
grangian system on the action Lie algebroid pr i :S\Gxg -> S\G. 

In addition, we may prove the following result 

Proposition 6 . 1 . (1) If •!> : TG —> S \ G x g is the map given by 

®(vh) = ( 7 r(/i), (T h l h - i(v h )), for all v h £ T h G (6.10) 

then $ is a fiberwise bijective Lie algebroid morphism over n. 

(2) The nonholonomic Lagrangian systems (L, D ) and (L, D) on TG and S/Gx 
0 are $ -related, that is, 

Io$ = I, 5(D) =5. 

(3) The system ( L , D) is regular and if 7 : I —> TG is a solution of the 
Lagrange-d Alembert equations for the system ( L,D ) then $ o 7 : I -t 
S\G x g is a solution of the Lagrange-dAlembert equations for the system 

(Z,D). 

Proof. (1) Consider the standard (right) action r of G on itself 

r : G x G G, ( h,h')£GxG -> r h> (h) = hh! £ G. 

As we know, the infinitesimal generator of r associated with an element to of g is 

log = Z 7 , 

where Xu is the left-invariant vector field on G such that ti7(e) = w. 

On the other hand, it is clear that the projection 7r : G —> S \ G is equivariant 
with respect to the actions r and ip. Thus, 

{T h n)(Xo{h)) = ^(co)(TT{h)), for h £ G. 
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Therefore, if p : S/G x g —> T{S\G) is the anchor map of the Lie algebroid 
pri : S\G x q —> S \ G, it follows that 

p($(tj(h))) = p(n(h),u) = (T h n)(uj(h)), for h G G. 

Furthermore, since 

< - 

[tu, u>'] = [u;,u/] 0 , for u>, u>' G 0 , 

we conclude that $ is a Lie algebroid morphism over 7 r. 

In addition, it is obvious that if h G G then 

5| T h G : T h G —► {7r(/l)} X 0 

is a linear isomorphism. 

(2) From (16.31) . (16.41) . (16.61) . (16.71) and (16.101) . we deduce that 

io<J> = L. 

Moreover, using (16.5[) . (| 6 . 8 j) . (16.91) and (I6.10L we obtain that 

$(£>) = D. 

(3) It follows from (1), (2) and using the results of Section [I] (see Theorem l4.6D . □ 

Next, we obtain the necessary and sufhcient conditions for a curve (h,u>) : I —> 
S' \ G x g to be a solution of the Lagrange-d’Alembert equations for the system 
(L, D). Let b< v > : 0 —> 0 * be the linear isomorphism induced by the scalar product 
<•,•>: 0 x 0 —> R and I : 0 —> 0 * be the inertia operator given by 

I(wi)(w 2 ) =< X(wi),w 2 >, for S 0 . ( 6 - 11 ) 

On the other hand, if b! G S \ G we will denote by : 0 —► T^,(S \ G) the linear 
epimorphism defined by 

\H^,(o/) = , T r (a 7 / )(/i / ), for u) 1 G 0 . 

In addition, if n(h') = h', we identify the vector space D-^, with the vector 
subspace D(h') of 0 . Then, using (13.71) . (16.71) and (|6.11l) . we deduce that the curve 
(h,u>) is a solution of the Lagrange-d’Alembert equations for the system (L,D) if 
and only if 

Kt) = ^ (t) (w(t)) 

(w(i) -5 1 (<<«%(*))) ^ 6 D± ty 

€ D h(ty 

for all t, where D~ is the orthogonal complement of the vector subspace C 0 
with respect to the scalar product < -, ■ >. These equations will be called the 
reduced Poincare- Chetaev equations. 

We treat next a simple example of the above general situation. 

The Veselova system. The most descriptive illustration of an LR system is the 
Veselova problem on the motion of a rigid body about a fixed point under the action 
of the nonholonomic constraint 

< ui, 7 >= 0 . 

Here, u> is the vector of the angular velocity in the body frame, 7 is a unit vector 
which is fixed in an space frame and < •, • > denotes the standard scalar product 
in R 3 (see jM] ). 

The Veselova system is an LR system on the Lie group G = SO(3) which is the 
configuration space of the rigid body motion. Thus, in this case, the Lie algebra 0 
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may be identified with R 3 and, under this identification, the Lie bracket [•, -] fl is the 
cross product x on R 3 . Moreover, the adjoint action of G = 50(3) on g = R 3 is the 
standard action of 50(3) on R 3 . This implies that < •, • > is an Adso( 3 )~invariant 
scalar product on g = R 3 . 

The vector subspace 0 of R 3 is just the orthogonal complement (with respect to 
< •, • >) of a vector subspace < 70 > of dimension 1, with 70 a unit vector in R , 
that is, 

5 = {w £ R 3 / < u, 70 >= 0 }. 

Therefore, 

s — 0^~ =< 70 > 

is a Lie subalgebra of g = R 3 . Furthermore, the isotropy group 5 of 70 with respect 
to the adjoint action of G = SO( 3), 

S = {s G SO(3)/s7<f = 7(f}, 

is a closed Lie subgroup with Lie algebra s. We remark that 5 is isomorphic to the 
circle 5 1 . 

Consequently, the corresponding homogeneous space M = S\SO( 3) is the orbit 
of the adjoint action of 50(3) on R 3 over the point 70 and, it is well-known that, 
such an orbit may be identified with the unit sphere 5 2 . In fact, the map 

5 \ 50(3) - 5 2 , [h] - 70 h = (h-^f 

is a diffeomorphism (see, for instance, EDI). 

Under the above identification the (right) action of SO(3) on M = S \ 50(3) is 
just the standard (right) action of 50(3) on 5 2 . Thus, our action Lie algebroid is 
the trivial vector bundle pr\ : 5 2 x R 3 —>■ 5 2 and the Lie algebroid structure on it 
is induced by the standard infinitesimal (right) action T : R 3 —> 3f(5 2 ) of the Lie 
algebra (R 3 , x) on 5 2 defined by 

’F(w)( 7 ) = 7 x u, for to g R 3 and 7 € 5 2 . 


In the presence of a potential V : 7 —► V (7) the nonholonomic Lagrangian system 
(L,D) on the Lie algebroid pry : 5 2 x R 3 —> R 3 is given by 

£(7.^) = ^I(w)(w) D( 7) = {7} x {w G R 3 / <u), 7 >= 0 }, 

I : R 3 —> R 3 being the inertia tensor of the rigid body. 

The Lagrange-d’Alembert equations for (L, D) are 


7 = 7 x w, w = I 1 {(Iw x ui) + u> x ——h A 7 }, < ui, 7 >= 0 (6-12) 

07 

where A is the Lagrange multiplier. Since the system (L, D) is regular, A is uniquely 
determined. In fact, 


A = — 


< lu) x w + 7 x fy,I _1 7 > 

< I - 1 7 ,7 > 


(6.13) 


Eqs (16.1211 and (16.131) are just the classical dynamical equations for the Veselova 
system (see [BB] ; see also [5B] IBB]). 
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6.3. Semidirect product symmetry and left action Lie algebroids. Here, we 
show how the reduction of some nonholonomic mechanical systems with semidirect 
product symmetry produces nonholonomic Lagrangian systems on left action Lie 
algebroids. 

Let us start by recalling the definition of a left action Lie algebroid (see [HU). Let 
(F, [•, -]f,Pf) be a Lie algebroid over a manifold N and tt : M —> N be a smooth 
map. A left action of F on tt : M —> N is a R-linear map 

'h : Sec(F) X{M), X G Sec(F) -> tf(X) G X(M) 

such that 

Mfx) = (/ ° *)*(*), *([*, y]f) = -MX), *00], 
(T m 7r)(^(X)(m)) = -p F (X(n(m))), 

for / G C°°(N), X, Y G Sec(F) and m G M. If * : Sec(E) —*■ X(M) is a left action 
of F on tt : M —> N and t f '■ F —> N is the vector bundle projection then the 
pullback vector bundle of F over tt, 

E = F*n = {(/, to) G F x M/r F (f) = i r(m)} 

is a Lie algebroid over M with Lie algebroid structure ([• > -}e ,Pe) which is charac¬ 
terized by 

[X, Y) e = [X, Y] f o tt, PE {X){m) = -$(X)(m), 

for X, Y G Sec(E) and m G M. The triple ( E, [■, -]b, p F ) is called the left action 
Lie algebroid of F over tt and it is denoted by F^tt (see [32]). 

Next, we consider a particular class of nonholonomic Lagrangian systems on 
left action Lie algebroids. Let V be a real vector space of finite dimension and 

• : G x V —> V be a left representation of a Lie group G on V. We also denote by 

• : g x V —> V the left infinitesimal representation of the Lie algebra g of G on V. 
Then, we can consider the semidirect Lie group S = G®V with the multiplication 

(g,v)(g',v') = (gg',v + g-v'). 

The Lie algebra s of S is the semidirect product s = gCDF with the Lie bracket 
s x s -> s given by 

\{u},v), {u',v')\ s = ([w,u/] 0 ,w -v' -J-v) 

for oj, to' G g and v, v' G V. Here, [•, -] 0 is the Lie bracket on g. 

Moreover, we use the following notation. If v G V then p v : g —> V is the linear 
map defined by 

p v {u) = oj ■ v, for w G g, 
and p* : V* —> g* is the dual map of p v : g —► V. 

Now, let N be a smooth manifold. Then, it is clear that the product manifold 
F = s x TN is the total space of a vector bundle over N. Moreover, if G s 

and X is a vector field on N then the pair ((ui,v),X) defines a section of the vector 
bundle t f : F = s x TN —> N. In fact, if {uij} is a basis of g, {uj} is a basis of 
V and {Xk} is a local basis of X(N) then {((oij, 0), 0), ((0, i)j), 0), ((0,0), Xfc)} is a 
local basis of Sec(F). 

The vector bundle t f : F —> N admits a Lie algebroid structure ([•, p F ), 
which is characterized by the following relations 

[((u,v),X),((u;',v'),X')} F = ([(cu,v),(cu',v%,[X,X'}) 

= ([w,w'j B ,w -v' -u' -v, [X,x']), 

p F {(u,v),X) = X, 

for ((w, *), X), ((a;', *'), X') G s x X(N). 


(6.14) 
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Next, suppose that i>o is a point of V and that O vo is the orbit of the action of 
G on V by vo, that is, 

O vo = {<?' G V/g G G}. 

Denote by n : M = N x O vo —► N the canonical projection on the first factor and 
by : Sec(F) — > X(M) the left action of F on 7 r, which is characterized by the 
following relation 

^((u i,u),X)(n,v) = (-X(n),u- v ) 

for ((uj,u),X) G s x X(N) and ( 71 , 17 ) G N x 0„ o = M. 

Then, we have the corresponding left action Lie algebroid t e ■ E = (s x 
TN)\;,tt —> M = TV x 0„ o . Note that E = (5 x TN)^n = (s x TW) x 0„ o and that 
the anchor map pE '■ E = (s x TN ) x 0„ o —> TM = TN x TO„ 0 of te : E —> M is 
given by 

p E {(w,u),X n ,v) = {X n , -oj ■ v) (6.15) 

for ((cj,u),X„,v) £sx T n N x O vo . 

Now, let L : (s x TN) x C>„ 0 — > R be a Lagrangian function and D be the vector 
subbundle of te '■ E —> M whose fiber over the point (n, v) G N x O vo = M 

is defined by 

D (n,v) = {(((w, u) ■ v),X n ), v)/u G g, X n G T n N } 

(6.16) 

C = (s x T n iV) x {u}. 

Next, we obtain the Lagrange-d’Alembert equations for the system (L, D). For this 
purpose, we choose a basis {ui a } of g, a basis {ua} of V, a system of local fibred 
coordinates (ad, id) on TN and a system of local coordinates (id) on O vo . Denote 
by l o a (respectively, u A ) the global coordinates on g (respectively, V) induced by 
the basis {w a } (respectively, {ma})- 
Suppose that 

[u7 a , g ^a/3^7 7 ^a ' 07 A — ^olA^B• 

Then, we have that 

r 1 n B - n C n B - n C n B 
C a/3 a iA ~ a pA a aC a aA a 0C- 

Next, we consider the local basis of sections {e*, e a , e^} of E given by 

ei(n,v) = ((0,0, ^r |n ), v), e a {n,v) = ((w a ,w a ■ v,0 n ),v) 

e A (n,v) = ((0,UA,0 n ),v) 


for (n,v) G N x O vo = M. Note that {e,,e Q } is a local basis of sections of the 
constraint subbundle D. In addition, if (ad, id ; y l , y a , y A ) are the local coordinates 
on E induced by the basis {e,, e a , e^}, it follows that 


V = x , y =uj 


V A = u A - cl a b UqUJ c 


(6.17) 


where u B is the local function on M = Nx O vo defined by u B = ( u B )\o V0 ■ Moreover, 


Ps(ei) = -g|i, PE(e a ) = Pa-^-r, p E (e A )=0, 

[Cq, , ep] E C a0 (e~, T d'yA^’O ^b): &a\e [^A, 7'Ck]e 07 a A^B , 


and the rest of the fundamental Lie brackets are zero. Thus, a curve 
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is a solution of the Lagrange-d’Alembert equations for the system ( L,D ) if and 
only if 

x 1 = y l , id = fd a y a , for all i and j. 
ddL dL 

iW 1 “3? = °' a “ *' 

d <1^> + + |4“)<&>»" - A S = 0, for all a, 


dt dy a 

y A = 0, for all A. 


If we consider the local expression of the curve in the coordinates (ad, id ; ad, w“, m a ) 
then, from (16.171) . we deduce that the above equations are equivalent to 

ad = y l , v J = for all * and j, 

d dL 

S<»> - aT = °' tor a “ 

d , dL . dL /o d , a p> dL . 

+ a^ c > + v(“ob"o 


dt 




<9L 


r)T 

+ 2a L U 0 d^A C l^ ~ P 1 ^ = °> fOT a11 


ii A = a A B UQLO a , for all A, 


or, in vector notation, 


v = —lo ■ v, 


ddL dL 
dt dx dx 
ddL ,*dL 

u = p v w. 


d dL *dL. * dL 

_ df ^ ” 2ad ^Pv-Q^) - Pvfo> 


Nonholonomic Lagrangian systems, of the above type, on the left action Lie alge- 
broid te : E = (s x TN) x O vo —> M = N x O vo may be obtained (by reduction) 
from an standard nonholonomic Lagrangian system with semidirect product sym¬ 
metry. 

In fact, let Q be the product manifold S x TV and suppose that we have a 
Lagrangian function L : TQ —* R and a distribution D on Q whose characteristic 
space ^((£,, v ),n) C T g GxT v V xT n N ~ T g GxVxT n N at the point ((<7, v), n) G SxN 
is 

D((g,v),n) = {((5)*))”) e TgG x V x T n N/v = (T g r g -i)(g) -v 0 }, (6.18) 

where vq is a fixed point of V. 

We can consider the natural left action of the Lie group 5 on Q and, thus, the left 
action A of the Lie subgroup H Vo = G Vo (S)V of S on Q , where G Vo is the isotropy 
group of Vo with respect to the action of G on V. The tangent lift TA of A is given 
by 

TA{{g,u), ( v g , (v,v),X n )) = {{T g l g )(v g ), (u + g-v,g- v),X n ) (6.19) 

for (g,u) G H Vo and (v g , (v,v),X„) G T(( g , v ),„)Q ^ T g G x V x T n N. 

Using (16.19|) , it follows that the distribution D is invariant under the action TA 
of H Vo on TQ. Moreover, we will assume that the Lagrangian function is also H Vo - 
invariant. Therefore, we have a nonholonomic Lagrangian system ( L,D ) on the 
standard Lie algebroid TQ —> Q which is H Vo -invariant. This type of systems were 
considered in [B3]. 
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Since the function L is 7?„ 0 -invariant, we deduce that there exists a real function 
L : (sxT N) x O vo —► R on the left action Lie algebroid te '■ E = (sxTN ) x 0„ o —► 
M = N x 0„ o which is defined by 

£(((w, *), ^n), w) = L((T e l g )(w), (v, g ■ v), X n ), (6.20) 

for (((w, v),X n ), u) £ s x T„7V x O vo , with g £ G, u = g~ 1 v$ and v £ V. 

Moreover, we may prove the following result. 

Proposition 6.2. (1) If $ : TQ ~ TG x (V x P) x TN -* £ = (s x TN) x O vo 

and p : G x V x N N x O vo are the maps defined by 

®{Ug,(V,V),X n ) = ((((Tglg-l^Ug )^- 1 ' i) , X n ) , g~ 1 ■ V Q ) , 

(b.zlj 

V(g,v,n) = (n, g -v 0 ), 

then & is a fiberwise bijective Lie algebroid morphism over ip. 

(2) The nonholonomic Lagrangian systems (L, D) and (L, D ) on TQ and E = 
(s x TN) x O vo are $ -related, that is, 

Lo$ = I, $(£))= D. 

Here, D is the vector subbundle of the vector bundle E whose fiber at the 
point (n,v) £ N x O vo is given by 116.16V . 

(3) If the system ( L,D ) is regular then the system ( L,D ) is also regular. In 
addition, if 7 : I TQ is a solution of the Lagrange-d Alembert equations 
for (L, D ) then $ o 7 : I —> (s x TN) x O vo is a solution of the Lagrange- 
d’Alembert equations for (L, D). 

Proof. (1) Suppose that u>\ and u >2 are elements of g, that v 1 and 62 are vectors of 
V and that X± and X 2 are vector fields on N. Then, we consider the vector fields 
Z\ and Z 2 on Q defined by 

Zi(g,v,n) = (u>i(g),g- Vi,X 1 (n)) £ T g G x V x T n N , 

Z 2 {g,v,n) = ( uj 2 (g),g ■ v 2 ,X 2 (n)) £ T g G x V x T n N , 

for ( g,v,n) £ G x V x N = Q, where Xj\ (respectively, Xv 2 ) is the left-invariant 
vector field on G such that tUi(e) = oj\ (respectively, tv 2 (e) = u) 2 ), e being the 
identity element of G. 

A direct computation proves that 

< - 

[Z 1 ,Z 2 )(g,v,n) = ([wi,W 2 ] 0 ( 5 ), 5 (wi -v 2 ~u 2 ■ t>i), [Xi,X 2 ](n)). 

Moreover, if ((wi,ui), X\) (respectively, ((u> 2 , v 2 ), X 2 )) is the section of the vector 
bundle te : E —> M induced by u> 1 , v 1 and X\ (respectively, u> 2 , v 2 and X 2 ) then 
it is clear that 

$ 0^1 = ((wi, ui), Xi) op, $oZ 2 = ({u 2 ,v 2 ),X 2 ) o p. 

Thus, using (16.141) . it follows that 

$ o [Zi, Z 2 \ = [((wi,*i),Xi), ({u> 2 ,i> 2 ), X 2 )]e o p. (6.22) 

On the other hand, we have that 

i. r L(^g V ^p){Ug,V,X n ) = ( X n , [Tgl g -l)(u g ) ■ (g ' VQ ) ) 

£ T n N x T g -i. Vo O Vo c T n N X V, 

for (g, v, n) £ Q and (u g ,v,X n ) £ T g G x V x T n N ~ T {g ^ v ^ n) Q. 

Therefore, from (16.151) and (16.211) . we deduce that 

Tp = pe 0 4>. 


(6.23) 
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Consequently, using (16.221) and (16.231) . we conclude that the pair (4>,yj) is a Lie 
algebroid morphism. Note that one may choose a local basis {Z*} of vector fields 
on Q such that 

Zi{g,v,n) = (t Ti(g),g- Vi,Xi(n)), for ( g,v,n ) € Q 

with Ui G 0, ii G V and Xi G X(N). 

Finally, if (g,v,n) G Q, it is clear that 

: T (g,v,n)Q ^ T g G x V x T n N -> E {n>g -i. Vo) ~0xFx T n N 
is a linear isomorphism. 

(2) It follows from (16.161) . (16.181) . (16.201) and (16.211) . 

(3) It follows using (1), (2) and the results of Section 0] (see Theorem 14.61) . □ 

The above theory may be applied to a particular example of a mechanical sys¬ 
tem: the Chaplygin Gyro (see [49[ [64]). This system consists of a Chaplygin 
sphere (that is, a ball with nonhomogeneous mass distribution) with a gyro-like 
mechanism, consisting of a gimbal and a pendulous mass, installed in it. The gim- 
bal is a circle-like structure such that its center coincides with the geometric center 
of the Chaplygin sphere. It is free to rotate about the axis connecting the north 
and south poles of the Chaplygin sphere. The pendulous mass can move along the 
smooth track of the gimbal. For this particular example, the vector space V is R 3 , 
the Lie group G is 50(3) and the manifold N is R 2 . The action of 50(3) on R 3 
is the standard one and vo = (0,0,1) is the advected parameter, see for more 
details. 

6.4. Chaplygin-type systems. A frequent situation is the following. Consider 
a constrained Lagrangian system ( L,D) on a Lie algebroid r: E —> M such that 
the restriction of the anchor to the constraint distribution, p\D: D —> TM, is an 
isomorphism of vector bundles. Let h: TM —> D C E be the right-inverse of p\ d, 
so that p o h = id tm- It follows that E is a transitive Lie algebroid and h is a 
splitting of the exact sequence 

0-^ Ker(p)- E —^ TM -^ 0 . 

Let us define the function L G C°°{TM) by L = L o h. The dynamics defined 
by L does not reduce to the dynamics defined by L because, while the map $ = p 
is a morphism of Lie algebroids and $(£)) = TM, we have Lo<b = Lohop^L. 
Nevertheless, we can use h to express the dynamics on TM, by finding relations 
between the dynamics defined by L and L. 

We need some auxiliary properties of the splitting h and its prolongation. We 
first notice that h is an admissible map over the identity in M, because pe 0 h = 
id tm and T idjvr °Ptm = *g?tm, but in general h is not a morphism. We can define 
the tensor K, a ker(p)-valued differential 2-form on M , by means of 

K(X, Y) = [ho X,hoY] — ho [X, Y] 

for every X, Y G X(M). It is easy to see that h is a morphism if and only if K = 0. 
In coordinates ( x l ) in M, ( x l ,v l ) in TM, and linear coordinates (x l ,y l ,y A ) on E 
corresponding to a local basis {e l ,e A } of sections of E adapted to the splitting h, 
we have that 

K = dx l A dx^ ® eA, 

2 J 

where are defined by [e*, ej] = 6 a- 

Since h is admissible, its prolongation T h h is a well-defined map from T(TM) 
to T e E. Moreover, it is an admissible map, which is a morphism if and only 
if ft, is a morphism. In what respect to the energy and the Cartan 1-form, we 
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have that ( T h h)*EL = E z and (T h h)*0L = 0 Z . Indeed, notice that by definition, 
( T h h)*EL = El oh and 

E L (h(v)) = ^ L(h(v) + t(h(v))\ t —o - L(h(v)) = ^ L(h(v + tv))\ t =o - L(h{v)) 

= ^ L(v + tv)\ t =o ~ L(v) = E l (v). 

On the other hand, for every V v = (v,w,V) £ T(TM ) = T™ (TM) where w = 
Tt(V), we have 

(( r h h)*e L , v) = (e L ,T h h(v, w, v)) = (e L , (h(v), h(w), th(v))) 

= |t=0 = ^j. L ( h ( V + tw )) |t=0 

= ^L(v + tw) |t=o = (6 l , V). 

Nevertheless, since h is not a morphism, and hence ( T h h )* o d ^ d o {[T h h )*, we 
have that ( T h h)*u)L ^ wy • Let JK be the 2-form on TM defined by 

JK v (V,W) = (J Kv) ,K Hv) (Tt m (V),Tt m (W))} 

where J is the momentum map defined by L and Kerp and V, W £ T h / V \(TM). 
The notation resembles the contraction of the momentum map J with the curvature 
tensor K. Instead of being symplectic, the map T h h satisfies 

( T h h)*u L =u L + JK. 

Indeed, we have that 

(T h h)*u> L -uj L = {do ( T h h )* - ( T h h)* o d] 0 L 

and on a pair of projectable vector fields U , V projecting onto X , Y respectively, 
one can easily prove that 

[d o C T h hf - (T h h)* O d) 0 L (U, V) = < e L , {T h h(U),T h h{V)\ - T h h([U , V])) 

from where the result follows by noticing that T h h o U is a projectable section and 
projects to hoX , and similarly T h hoV projects to hoY. Hence [T h h(U), T h h(V )] — 
T h h([U, V] is projectable and projects to K (X, Y). 

Let now T be the solution of the nonholonomic dynamics for ( L 7 D) 7 so that L 
satisfies the equation irWL — dE l £ D° and the tangency condition r| D £ T D D. 
From this second condition we deduce the existence of a vector field f £ X(TM) 
such that T h hoT = Toh. Explicitly, the vector field f is defined by F = T p poToh , 
from where it immediately follows that F is a SODE vector field on M. 

Taking the pullback by T h h of the first equation we get ( T h h )* {iv^L — dE^) = 0 
since ( T h h)*D° = 0. Therefore 

0 = ( T h h)*i r uj L - (T h h)*dE L 
= i T (T h h)*w L - d(T h h)*E L 
= if (u>l + JK) — dEi 
= Zf — dEi ifJK. 

Therefore, the vector field f is determined by the equations 

^ru JZ -dE L = -(J,K(T, ■)), 

where T is the identity in TM considered as a vector field along the tangent bundle 
projection tm (also known as the total time derivative operator). Equivalently we 
can write these equations in the form 

drd L -dl=(J,K(T 1 ■)). 
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Note that if a : / -t TM is an integral curve of f then a = h o a : I —> D is 
a solution of the constrained dynamics for the nonholonomic Lagrangian system 
(L, D ) on E. Conversely, if a : I —► D is a solution of the constrained dynamics 
then po a : I -> TM is an integral curve of the vector field f. 

Finally we mention that extension of the above decomposition for non transitive 
Lie algebroids is under development. 

Chaplygin systems and Atiyah algebroids. A particular case of the above 
theory is that of ordinary Chaplygin systems(see [U [7] 13 [37] and references there 
in). In such case we have a principal G-bundle n: Q —* M = Q/G. Then, we 
may consider the quotient vector bundle E = TQ/G —> M = Q/G and, it is 
well-known that, the space of sections of this vector bundle may be identified with 
the set of G- invariant vector fields on Q. Thus, using that the Lie bracket of two 
G-invariant vector fields is also G-invariant and the fact that a G-invariant vector 
field is 7r-projectable, we may define a Lie algebroid structure ([•,•], p) on the vector 
bundle E = TQ/G —> M = Q/G. The resultant Lie algebroid is called the Atiyah 
(gauge) algebroid associated with the principal bundle tt : Q —► M = Q/G (see 
|48| ). Note that the canonical projection d* : TQ —► E = TQ/G is a fiberwise 
bijective Lie algebroid morphism. Now, suppose that ( Lq,Dq ) is an standard 
nonholonomic Lagrangian system on TQ such that Lq is G-invariant and Dq is the 
horizontal distribution of a principal connection on 7r : Q — > M = Q/G. Then, we 
have a reduced nonholonomic Lagrangian system (L, D ) on E. In fact, Lq = L o <f> 
and d>((H(g) 9 ) = -D^q), for all q £ Q. Moreover, p\ D : D —+ TM = T(Q/G) is 
an isomorphism (over the identity of M) between the vector bundles D —> M and 
TM -» M. Therefore, we may apply the above general theory. 

Next, we describe the nonholonomic Lagrangian system on the Atiyah algebroid 
associated with a particular example of a Chaplygin system: a two-wheeled planar 
mobile robot (see m and the references there in). Consider the motion of two¬ 
wheeled planar mobile robot which is able to move in the direction in which it points 
and, in addition, can spin about a vertical axis. Let P be the intersection point of 
the horizontal symmetry axis of the robot and the horizontal line connecting the 
centers of the two wheels. The position and orientation of the robot is determined, 
with respect to a fixed Cartesian reference frame by (x, y, 9) £ SE( 2), where 6 £ S 1 
is the heading angle, the coordinates (x,y) £ M 2 locate the point P and SE( 2) is 
the group of Euclidean motions of the two-dimensional plane R 2 . Let 'i/i , £ 

S 1 denote the rotation angles of the wheels which are assumed to be controlled 
independently and roll without slipping on the floor. The configuration space of 
the system is Q = T 2 x SE( 2), where T 2 is the real torus of dimension 2. 

The Lagrangian function Lq is the kinetic energy corresponding to the metric 
9Q 

gQ = mdx ® dx + mdy <8> dy + mol cos 9{dy ® dO + dO ®> dy ) 

—mol sin 9{dx <g> dO + dO ® dx) + Jdd ® dO + ^dt/i ® di/i + Jid^i ® di/> 2, 

where m = mo + 2m\, mo is the mass of the robot without the wheels, J its 
momenta of inertia with respect to the vertical axis, m\ the mass of each wheel, J 2 
the axial moments of inertia of the wheels, and l the distance between the center 
of mass G of the robot and the point P. Thus, 

Lq = — (mx 2 + my 2 + 2moly9 cos 9 — 2molx9 sin 9 

+J9 2 + J 21 P 2 + hipl)- 

The constraints, induced by the conditions that there is no lateral sliding of the 
robot and that the motion of the wheels also consists of a rolling without sliding, 
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are 


x sin 9 — y cos 9 = 0, 
x cos 9 + y sin 9 + c9 + Ripi = 0, 
xcos9 + y sin9 — c9 + Rip 2 = 0, 


where R is the radius of the wheels and 2c the lateral length of the robot. The 
constraint distribution D is then spanned by 


{Hi 


d 

Wi 


R f a d 

2 (C0S ^ 


sin0 


d Id 


), 


dy cd9' 

„ 9 R < H 9 


. ,, d 13 


Note that if {6,6, £ 3 } is the canonical basis of se(2), 


[£1.62] =0, [6,6] =-6> [6,6] = 6, 


then 


Hi 


d 
dip 1 




H 2 


d R<- R<- 

d^- 1 2^ + 2c 


where 6 (i = 1, 2, 3) is the left-invariant vector field of SE( 2) such that 6 (e) = 6, 
e being the identity element of SE( 2). 

On the other hand, it is clear that Q = T 2 x SE( 2) is the total space of a 
trivial principal 5£'(2)-bundle over M = T 2 . Moreover, the metric gQ is SE( 2)- 
invariant and Dq is the horizontal distribution of a principal connection on Q = 
T 2 x SE(2) -► T 2 . 

Now, we consider the corresponding Atiyah algebroid 


E = TQ/SE{ 2) ~ (IT 2 x TSE(2))/SE(2) —> M = T 2 . 


Using the left-translations on SE( 2), we have that the tangent bundle to SE(2) may 
be identified with the product manifold SE( 2) xse(2) and, under this identification, 
the Atiyah algebroid is isomorphic to the trivial vector bundle 

f T 2 = r T 2 opn : TT 2 x se(2) —► T 2 , 


where r T 2 : TT 2 —► T 2 is the canonical projection. In addition, if ([•,•], p) is the 

d d 

Lie algebroid structure on t T 2 : TT 2 x se(2) —► T 2 and { 77 —, =— , 6 , 6 , 6 } is the 

dip 1 dip 2 

canonical basis of sections of f T 2 : TT 2 x se(2) —> T 2 then 


d 


"W = 


P(6) = 0, * = 1,2,3 


d <9 d 

Wi P ^ = dih’ 

[ 6 , 6 ] = - 6 , [ 6 , 6 ] = 6 , 

and the rest of the fundamental Lie brackets are zero. 

Denote by (V’i,V’ 2 ,Vh,w 1 ,w 2 ,w 3 ) the (local) coordinates on TT 2 x se(2) 
d d 

induced by the basis {vrt - , vrt - , 6 , 6 , 61 - Then, the reduced Lagrangian L : 
dipi dip 2 

TT 2 x se(2) —► R is given by 


L = -( m (w 1 ) 2 + m(w 2 ) 2 + 2itioIuj 2 uj 3 + J(w 3 ) 2 + 6*/) 2 + J 262 ) 


and the constraint vector subbundle D is generated by the sections 


d 

ei ~Wi 




e-2 


d R R 
d^ ~ 2^ + 2c^ 3 ' 


Since the system ( Lq,Dq ) is regular on the standard Lie algebroid tq : TQ —> Q , 
we deduce that the nonholonomic Lagrangian system (L, D) on the Atiyah algebroid 
Tf 2 : TT 2 x se(2) —► T 2 is also regular. 
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Now, as in Section [3j we consider a basis of sections of f T 2 : T T 2 x se(2) —► T 2 
which is adapted to the constraint subbundle D. This basis is 

( e l> e 2, £ 1 , £ 2 , £ 3 }- 

The corresponding (local) coordinates on T T 2 xse(2) are (^ 1 , 4>2, V’i, Va w 1 , w 2 , <I> 3 ), 
where 

1 ~ 1 P : P J 2 ~ 2 3 ~ 3 P : P : 

U> =UJ - —ipi - —1P2, U) =OJ , U> =U - —4> 1 + —1P2- 
2 2 2c 2c 

Therefore, using (13TD . we deduce that the Lagrange-d’Alembert equations for the 

system ( L,D ) are 

V U(lp2-1pl) rr> : , C ! x 7 U(lp 2 -1pl) , , 

V’i = p‘2 _ S 2 (Pfa+Slp 1 ), V’2 =- p 2 _ 5-2 ( P V>1 + SV’2), 

w 1 = w 2 = w 3 = 0, 
where P, S' and [/ are the real numbers 

„ P 2 . J, T 0 P 2 , J, Tr R 3 , 

P = -( to + -) + j 2 , 5 = _( to __) ; U = — moL 
On the other hand, the Lagrangian function L : TT 2 —> K on TT 2 is given by 

L{ip 1 ,^ 2 , V’i, V 2 ) = + PV 2 + 2SV’iV’2) 

and the 1-form (J, PT(T, •)) on TT 2 is 

(J,PT(T, ■)) = -U(ip2 - V’i)(V’i#2 - V>2<# 1 ). 


7. Nonlinearly constrained Lagrangian systems 


We show in this section how the main results for linearly constrained Lagrangian 
systems can be extended to the case of Lagrangian systems with nonlinear nonholo- 
nomic constraints. This is true under the assumption that a suitable version of the 
classical Chetaev’s principle in nonholonomic mechanics is valid (see e.g., jjl] for 
the study of standard nonholonomic Lagrangian systems subject to nonlinear con¬ 
straints). 

Let r : E —► M be a Lie algebroid and A4 be a submanifold of E such that 
tt = t : At —> M is a fibration. At is the constraint submanifold. Since n is a 
fibration, the prolongation T E A4 is well-defined. We will denote by r the dimension 
of the fibers of n: At —► M, that is, r = dim At — dim M. 

We define the bundle V —> A4 of virtual displacements as the subbundle of 
t*E of rank r whose fiber at a point a £ At is 


V a = {b£ E T{a) | b v a G T a M }. 


In other words, the elements of V are pairs of elements (a, b) G E © E such that 


dt 


(j>{a + tb) 


t =0 


= 0, 


for every local constraint function (j). 

We also define the bundle of constraint forces \& by U/ = S*((T E M)°), in 
terms of which we set the Lagrange-d’Alcmbert equations for a regular Lagrangian 
function L G C°°(E) as follows: 


(j r wi - dE L )\ M G Sec(’J'), 
T| M G Sec(T E JA), 


the unknown being the section T. The above equations reproduce the corresponding 
ones for standard nonlinear constrained systems. 
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From ( 12 . 81 ) and ( 17 . 11 ) . it follows that 

(isrUL - *awl)|a4 = 0, 

which implies that a solution T of equations (EH) is a SODE section along Ad, that 
is, (ST — A)|Ad = 0. 

Note that the rank of the vector bundle (T E Ad)° —> Ad is s = rank (S') — r and, 
since n is a fibration, the transformation S* : (T E JA)° —► T defines an isomorphism 
between the vector bundles (T E A4)° —> Ad and T —> Ad. Therefore, the rank of T 
is also s. Moreover, if a £ Ai we have 

= S*((T E Ad)°) = { C o Tr | C £ V a ° }. (7.2) 

In fact, if a a £ (T E Ad)°, we may define ( £ ^r(a) by 

((b) = a a (( v (a, 6)), for b £ E r(a) . 

Then, a direct computation proves that ( £ V° and S*(a a ) = ( oTt. Thus, we 
obtain 

* a C{(oTT\(eV° a } 

and, using that the dimension of both spaces is s, we deduce ( 17 . 21 ) holds. Note 
that, in the particular case when the constraints are linear, we have V = t*(D) and 
T = D°. 

Next, we consider the vector bundles F and T V M over A4 whose fibers at the 
point a £ M are 

F a = cj^i^a), r^M = { (b,v) £ V a x T a M \ Tn(v) = p(b) }. 

It follows that 

F a = [z £ T E E | exists (eV° and i z UL(a) = ( o Tt } 

and 

T a v M = { * £ T e M | Tn(z) £ V Q } = { z £ T E M \ S(z) £ T E M }. (7.3) 

Note that the dimension of M is 2r and, when the constraints are linear, i.e., 
AI is a vector subbundle D of E, we obtain 

T^M = T a °D , for all a £ M = D. 

Moreover, from (17731) . we deduce that the vertical lift of an element of V is an 
element of T v Ad. Thus we can define for b, c £ V a 

Ga’ v (b , c) = uj L (a)(b, ( v (a, c)), 

where b £ T e E and Tr(b) = b. 


7.1. Dynamics in local coordinates. Here we analyze the local nature of equa¬ 
tions EH- We consider local coordinates ( x l ) on an open subset U of M and take 
a basis {e a } of local sections of E. In this way, we have local coordinates (ad,y a ) 
on E. Suppose that the local equations defining Ad as a submanifold of E are 

<t> A = 0, A = 1,..., s, 


where < f> A are independent local constraint functions. Since tt : A4 —> M is & 

dd>A 

fibration, it follows that the matrix (——) is of rank s. Thus, if d is the differential 

K dy a) 

of the Lie algebroid T e E —> E , we deduce that {g?</>' 4 |a 4 }a=i,...,s is a local basis of 
sections of the vector bundle (T E Ad)° —> Ad. Note that 
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QlA 

Moreover, {S*{dcj) A )\M = 77 — X a \M}A=i,...,s is a local basis of sections of the 


dy 

vector bundle 4/ —> At. 

Next, we introduce the local sections {Za}a= i,. 


of T e E -a E defined by 


%z a ul = S*{d<f A ) = 

A direct computation, using (ESI), proves that 

d(f A 


~dy° 


■X° 


Z A = - 


dy c 


-W a0 Vg 


for all A, 


d 2 L 


(7.4) 


where (W a0 ) is the inverse matrix of (W a g = q a. fs )• Furthermore, it is clear 

that {Za\m} is a local basis of sections of the vector bundle F —> At. 

On the other hand, if T/, is the Euler-Lagrange section associated with the regular 
Lagrangian L, then a section T of T E At —> A4 is a solution of equations ED if 
and only if 

r = (t l + x a z a )\m 

with A " 4 local real functions on E satisfying 

(. \ A dcj ) B {ZA)+d ( t > B {T L ))\M = 0 , for all B = 1 ,..., s. 

Therefore, using ED, we conclude that there exists a unique solution of the 
Lagrange-d’Alembert equations ED if and only if the matrix 


(C AB = 


d<j/ 


W 


oi(3 


d(j> 1 


') A.B= 1 . 


dy a '' dyP 

is regular. We are now ready to prove the following result. 


(7.5) 


Theorem 7.1. The following properties are equivalent: 

(1) The constrained Lagrangian system (L,A4) is regular, that is, there exists 
a unique solution of the Lagrange-d’Alembert equations, 

(2) KerG L ' v = {0}, 

(3) T E MnF= {0}, 

(4) T v Mn(T v M) E = {0}. 

Proof. It is clear that the matrix ( C AB ) in (17.51) is regular if and only if T E A4nF = 
{0}. Thus, the properties (1) and (3) are equivalent. Moreover, proceeding as in 
the proof of Theorem 13. 121 we deduce that the properties (2) and (3) (respectively, 
(2) and (4)) also are equivalent. □ 

Remark 7.2 (Lagrangians of mechanical type). If L is a Lagrangian function of 
mechanical type, then, using Theorem 17.11 we deduce (as in the case of linear 
constraints) that the constrained system (L , M) is always regular. • 

7.2. Lagrange-d’Alembert solutions and nonholonomic bracket. Assume 
that the constrained Lagrangian system (L, At) is regular. Then (3) in Theorem 1 7. II 
is equivalent to [T e E)\m = T E A4 ® F. Denote by P and Q the complementary 
projectors defined by this decomposition 

P a :T E E^T E M, Q a :T E E-^F a , for all a G M. 

As in the case of linear constraints, we may prove the following. 

Theorem 7.3. Let (L, Ai) be a regular constrained Lagrangian system and letT l 
be the solution of the free dynamics, i.e., ir L ajL = dE^. Then, the solution of the 
constrained dynamics is the SODE T(l,m) obtained as follows 

F (l,m) = P(Fl\m)- 
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On the other hand, (4) in Theorem 17.11 is equivalent to {T e E)\j V [ = T V A4 © 
(T v M) 1 - and we will denote by P and Q the corresponding projectors induced by 
this decomposition, that is, 

P a : T 0 e E -► r a v A4, Q a : T e E -> (T a v M) x , for all a G M. 

Theorem 7.4. Let (L, JVf) be a regular constrained Lagrangian system, Tl (respec¬ 
tively, be the solution of the free (respectively, constrained) dynamics and 

A be the Liouville section of T e E —> E. Then, T(l,m) = P(Tl\m) If and only if 
the restriction to A4 of the vector field p 1 (A) on E is tangent to A4. 

Proof. Proceeding as in the proof of Lemma 13.111 we obtain that 

0 Ver (T a E E) = F a , for all a G M. 

Thus, it is clear that 

Q(T l (o)) G F a C for all a G M. 

Moreover, from (17.31) and using the fact that the solution of the constrained 
namics is a SODE along M, we deduce 

r (L,A 4 )(a) = P(^L(a)) G T a v Al, for all a G M, 

if and only if the restriction to A4 of the vector field p 1 (A) on E is tangent to 
This proves the result. 


dy- 

M. 

□ 


Remark 7.5 (Linear constraints). Note that if A1 is a vector subbundle D of E, 
then the vector field p 1 (A) is always tangent to A i = D. • 

As in the case of linear constraints, one may develop the distributional approach 
in order to obtain the solution of the constrained dynamics. In fact, if (L,M) is 
regular, then T v M —► M is a symplectic subbundle of ( T e E,ujl ) and, thus, the 
restriction uj L:M of lul to T V A4 is a symplectic section on that bundle. We may 
also define s L ’ M as the restriction of dE l to T v A4. Then, taking the restriction of 
Lagrange-d’Alembert equations to T v M, we get the following equation 

i r cv L ’ M = e L ’ M , (7.6) 

which uniquely determines a section T of T V A4 —*■ M. It is not difficult to prove 
that f = P{Tl\m)- Thus, the unique solution of equation ( 17 . 61 ) is the solution of 
the constrained dynamics if and only if the vector field p 1 (A) is tangent to M. 

Let (L,M) a regular constrained Lagrangian system. Since S* : ( T E M)° — > T 
is a vector bundle isomorphism, it follows that there exists a unique section j^) 
of ( T e M)° -> M such that 

i Q(T L \ M )OJL = S'* (<*(!,>1))- 
Moreover, we have the following result. 


Theorem 7.6. If (L, M) is a regular constrained Lagrangian system and T^ is 
the solution of the dynamics, then dr (i m) (El\m) = 0 if and only if cx(l,m)(^\m) = 
0. In particular, if the vector field p 1 (A) is tangent to A4 , then dr (L m) (El\m) = 0- 

Proof. From Theorem 17.31 we deduce 


0 i r iL ,M) u) L - dE L )\ M = —S*(a( LtM )). 
Therefore, using that T (l,m) is a SODE along M, we obtain 

d r (L ,M)( E L\M) = U( L , M )(A\m)- 


□ 
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Now, let (L,Ad) be a regular constrained Lagrangian system. In addition, sup¬ 
pose that / and g are two smooth functions on Ad and take arbitrary extensions 
to E denoted by the same letters. Then, as in Section 13.51 we may define the 
nonholonomic bracket of / and g as follows 

{f,g}nh = UL(P(X f ),P(X g ))\ M , 

where Xf and X g are the Hamiltonian sections on T e E associated with / and g, 
respectively. 

Moreover, proceeding as in the case of linear constraints, one can prove that 

f = p\R L )(f) + {f,E L } nh , f G C°°(M), 

where Rl is the section of T E M —> Ad defined by Rl = P{Tl\m) ~ P{^l\m)- 
Thus, in the particular case when the restriction to Ad of the vector field p 1 (A) on 
E is tangent to .Ad, it follows that 

f={f,E L } nh , for feC°°(M). 

Alternatively, since T v Ad is an anchored vector bundle, we may consider the 
differential df G Sec((T v Ad)*) for a function / G Thus, since the restric¬ 

tion lo L ’ M of wl to T v Ad is regular, we have a unique section Xf G Sec(T v Ad) 
given by ix f u> L,M = df and it follows that 

{f,g}nh = u LM (X f ,X g )- 

7.3. Morphisms and reduction. Let ( L , Ad) be a regular constrained Lagrangian 
system on a Lie algebroid r : E —> M and let (L',Ad ') be another constrained 
Lagrangian system on a second Lie algebroid t' : E' —>■ M’. Suppose also that we 
have a fiberwise surjective morphism of Lie algebroids d* : E — > E' over a surjective 
submersion cf> : M —>■ M' such that: 

(i) L = L' o$, 

(ii) $|a 4 : M M' is a surjective submersion, 

(hi) $(V a ) = V4 (a) , for all a G M. 

Note that condition (ii) implies that d>(V a ) C for all a G Ad. Moreover, if 

V ($) is the vertical bundle of d> and 

V r Q ($) C T a M, for all a G M, 

then condition (ii) also implies that V4( a ) C $(V a ), for all a G Ad. 

On the other hand, using condition (iii) and Proposition 14.21 it follows that 
kerG L ,v = {0} and, thus, the constrained Lagrangian system (L', Ad') is regular. 
Moreover, proceeding as in the proof of Lemma 14.51 and Theorem 14.61 we deduce 
the following results. 

Lemma 7.7. With respect to the decompositions 

(■ T e E)\ m = R E M © F and (T E 'E')\ M , = T E> M' 0 F' 
we have the following properties 

( 1 ) =T E 'M', 

(2) T $ $(F) = F', 

(3) If P, Q and P\Q’ are the projectors associated with (L, Ad) and (L',Ad'), 

respectively, then P' o T $ <f> = o P and Q' o o Q. 

With respect to the decompositions 

(■ T e E)\ m = T V M 0 (T v M) e and {T E ’E')\ M . = T v M' © (T v M') E 
we have the following properties 

(4) (T $ d>)(T v Ad) = T v 'Ad', 
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(5) =_ (T V M')\ 

(6) If P, Q and P' and Q' are the projectors associated with ( L , M) and (L 1 , M'), 
respectively, then P' o = T $< f> o P and Qf o T*<I> = o Q. 

Theorem 7.8 (Reduction of the constrained dynamics). Let (L,M) be a regular 
constrained Lagrangian system on a Lie algebroid E and let (L 1 , Ad') be a con¬ 
strained Lagrangian system on a second Lie algebroid E'. Assume that we have 
a fiberwise surjective morphism of Lie algebroids <£> : E —> E' over (f> : M —> M' 
such that conditions (i)-(iii) hold. IfU^ LM ^ is the constrained dynamics for L 
andT (£/ is the constrained dynamics for L', respectively, then ^ = 

r(L',A 4 ') 0 If aft) is a solution of Lagrange-d’Alembert differential equations for 
L, then <f>(a(f)) is a solution of Lagrange-d’Alembert differential equations for L'. 

We will say that the constrained dynamics Iis the reduction of the 
constrained dynamics L(j^ M , by the morphism $. As in the case of linear 
constraints (see Theorem 14. 71) . we also may prove the following result 

Theorem 7.9. Under the same hypotheses as in Theorem | 7. 6] we have that 
if o<S>,g'o<I>} nh = {f',g'y nh o<f>, 

for f',g' e C°°(M'), where {•, •} n h (respectively, is the nonholonomic 

bracket for the constrained system (L,M) (respectively, (L',M')). In other words, 
d> : Ad —> M' is an almost-Poisson morphism. 

Now, let </> : Q —> M be a principal G-bundle and r : E — > Q be a Lie algebroid 
over Q. In addition, assume that we have an action of G on E such that the quotient 
vector bundle E/G is defined and the set Sec(-E)° of equivariant sections of I? is a 
Lie subalgebra of Sec (E). Then, E' = E/G has a canonical Lie algebroid structure 
over M such that the canonical projection $ : E —> E' is a fiberwise bijective Lie 
algebroid morphism over </> (see Theorem 14.811 . 

Next, suppose that (L,M) is a G-invariant regular constrained Lagrangian sys¬ 
tem, that is, the Lagrangian function L and the constraint submanifold A4 are 
G-invariant. Then, one may define a Lagrangian function L' : E' —> R on E' such 
that 

L = L' o$. 

Moreover, G acts on A4 and if the set of orbits A4' = M./G of this action is a 
quotient manifold, that is, MI is a smooth manifold and the canonical projection 
d ) |m '■ M = M/G is a submersion, then one may consider the constrained 

Lagrangian system {L',M') on E'. 

Remark 7.10 (Quotient manifold). If M is a closed submanifold of E, then, 
using a well-known result (see |ll Theorem 4.1.20]), it follows that the set of orbits 
MI = M/G is a quotient manifold. • 

Since the orbits of the action of G on E are the fibers of d> and M is G-invariant, 
we deduce that 

V a ($) C T a M, for all a G M, 

which implies that $|v a : V a —» V^( a ) is a linear isomorphism, for all a £ M. 

Thus, from Theorem El we conclude that the constrained Lagrangian system 
( L',M' ) is regular and that 

T $ <i) or (IiM) = r° $, 

where ^ (resp., ^ q) is the constrained dynamics for L (resp., L'). In 
addition, using Theorem 17.91 we obtain that $ : M —> MI is an almost-Poisson 
morphism when on M and M' we consider the almost-Poisson structures induced 
by the corresponding nonholonomic brackets. 

We illustrate the results above in a particular example in the following subsection. 
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7.4. Example: a ball rolling on a rotating table. The following example is 
taken from [U HU 37 ] . A (homogeneous) sphere of radius r > 0, unit mass m = 1 
and inertia about any axis k 2 , rolls without sliding on a horizontal table which 
rotates with constant angular velocity fl about a vertical axis through one of its 
points. Apart from the constant gravitational force, no other external forces are 
assumed to act on the sphere. 

Choose a Cartesian reference frame with origin at the center of rotation of the 
table and z-axis along the rotation axis. Let (x, y ) denote the position of the point 
of contact of the sphere with the table. The configuration space for the sphere on 
the table is Q = R 2 x S'O(S), where SO( 3) may be parameterized by the Eulerian 
angles 9, ip and ip. The kinetic energy of the sphere is then given by 

T = i(i; 2 + y 2 + k 2 (9 2 + ip 2 + 2ipipcos9)), 

and with the potential energy being constant, we may put V = 0. The constraint 
equations are 

x — r9smip + rip sin 9 cos ip = —fly, 
y + r9 cos ip + rip sin 9 sin ip = fhc. 

Since the Lagrangian function is of mechanical type, the constrained system is 
regular. Note that the constraints are affine, and hence not linear, and that the 
restriction to the constraint submanifold Ai of the Liouville vector field on TQ is 
not tangent to A4. Indeed, the constraints are linear if and only if SI = 0. 

Now, we can proceed from here to construct to equations of motion of the sphere, 
following the general theory. However, the use of the Eulerian angles as part of the 
coordinates leads to very complicated expressions. Instead, one may choose to 
exploit the symmetry of the problem, and one way to do this is by the use of 
appropriate quasi-coordinates (see m- First of all, observe that the kinetic 
energy may be expressed as 

T = -(x 2 +y 2 + k 2 {io 2 + to 2 + w 2 )), 

where 

uj x = 9 cos ip + ip sin 9 sin ip, 
bjy = 9 sin ip — ip sin 9 cos ip, 
u) z = ip cos 9 + ip, 

are the components of the angular velocity of the sphere. The constraint equations 
expressing the rolling conditions can be rewritten as 

x — rw y = —fly, 
y + ru> x = f lx. 

Next, following m, we consider local coordinates (x, y, 9, ip, ip-, 7r*)* = i > „. J 5 on TQ = 
TIR 2 x T(SO( 3)), where 

x = x, y = y, 9 = 9, (p = ip, ip = ip, 

7Ti = rx + k 2 q 2 , tt 2 = ry - k 2 qi, -k 3 = k 2 q 3 , 

k 2 k 2 

^4= ( fc2+r 2 j (±-rq 2 + Sly), tt 5 = (y + rh-Slx), 

and ( 51 , 92 , 93 ) are the quasi-coordinates defined by 

91 = u x , 92 = u y , 93 = w z . 
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As is well-known, the coordinates qi only have a symbolic meaning. In fact, 
d d d 

{ ——, ——, ——} is the basis of left-invariant vector fields on SO( 3) given by 
oqi dq 2 dq 3 


8 i n 3 
^- = (cos«- . sin# , av 


sin 'll) , d n d . 


_d_ 

dqi 


d . d cos ip d d 

d d 

dq 3 dip ’ 

and we have that 

<9 9 d d <9 d d d 

Note that in the new coordinates the local equations defining the constraint 
submanifold M are 7 T 4 = 0, 755 = 0. On the other hand, if P : (T t< ^TQ)\m = 
Tm(TQ) —> T t ®M = TM and Q : Tm{TQ) —> F are the projectors associated 
with the decomposition Tji 4 (TQ) = TM. ® F, then we have that (see [TT] ) 

d <9 

Q = w — <8 a7r4 + -— ® dn 5 , 

07T4 C/TTs 

d d 

P = Id — —— ® d 7 r 4 - —— (g) dn 5 . 

OIT 4 (77T5 

Moreover, using that the unconstrained dynamics is given by 

. d . d ■ d . d ■ d 

L X dx y dy dd ^dtp dip (fc 2 + r 2 ) y i 97T4 ' (fc 2 + r 2 ) c?7r 5 

. <9 . d . d . d . d fc 2 0 . d fc 2 0 . d 

91 <9<7i + dq 2 ^ ^ dq 3 + (fc 2 + r 2 ) ^£>774 (k 2 + r 2 ) X dn 5 ' 


k 2 fl d 
:V- 


k 2 fl d 

-x- 


= X d^ +V d-y 


we deduce that the constrained dynamics is the SODE F( L along M defined by 


= (Pril "> = (± S + 8 Yy + * S + + ^ )U 


.. d . d . d . d . <9 

= (itt + 2/wr +91-5-f 92-h 93-^—) At- 

dx dy dqi dq 2 dq 3 


(7.7) 


This implies that 


d r (L , M) (E L \m) = dr, 


Q 2 k 2 

(l, M ) ( l \m ) = Tprr-ji~( x± + yy)U- 


(fc 2 + r 2 ) 

Consequently, the Lagrangian energy is a constant of the motion if and only if 
0 = 0 . 

When constructing the nonholonomic bracket on A4, we find that the only non¬ 
zero fundamental brackets are 

{x, 7Ti }nh = r, { y , 7T2 }nh = T j 

{9l) = —1, {92, Tl}n/i = 1, {93: 7T 3 }n7i = 1, 

fc 2 rfc 2 0 /y ol 

{Tl, 7T2}n/i — 77"3) I 71 " 2 ? ^ jn/i — ^2)^2 _|_ r 2^’ ' ' 

k 2 rk 2 fl 

{7T 3 ; 7Tl}ra/i = (fc 2 +r 2^ 2 “ (p+r 2 )*’ 

in which the “appropriate operational” meaning has to be attached to the quasi¬ 
coordinates qi . As a result, we have 

f = R L (f) + {f,L} nh , for / G C°°(M) 
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where Rl is the vector field on A i given by 


Rl 


k 2 n d d rfl d d 

( k 2 +r 2 ) X dy ^dx (k 2 + r 2 ) X dq\ + ^ dq 2 

d d d 

+ x(n 3 - k 2 fl)- -h £/(tt 3 - k 2 fl)- - k 2 (TTiX + -K 2 y)^—))\M- 

OTTl O 7T2 07T3 


Note that Rl = 0 if and only if Q = 0. 

Now, it is clear that Q = R 2 x SO( 3) is the total space of a trivial principal 50(3)- 
bundle over R 2 and the bundle projection (/> : Q —> M = R 2 is just the canonical 
projection on the first factor. Therefore, we may consider the corresponding Atiyah 
algebroid E' = TQ/SO( 3) over M = R 2 . Next, we describe this Lie algebroid. 

Using the left-translations in 50(3), one may define a diffeomorphism A between 
the tangent bundle to 50(3) and the product manifold 50(3) x R 3 (see pQ). In 
fact, in terms of the Euler angles, the diffeomorphism A is given by 


A(6>, ip, ip; 0, p, ip) = (0, p, ip; u; x ,ujy,v z ). (7.9) 

Under this identification between T(SO( 3)) and 50(3) x R 3 , the tangent action of 
50(3) on T(50(3)) = 50(3) x R 3 is the trivial action 

50(3) x (50(3) x R 3 ) —> 50(3) x R 3 , (g, (h,u)) ha (gh, w). (7.10) 


Thus, the Atiyah algebroid TQ/SO( 3) is isomorphic to the product manifold TR 2 x 
R 3 , and the vector bundle projection is 7 r2 opri, where pr\ : TR 2 x R 3 —> TR 2 and 
t R 2 : TR 2 R 2 are the canonical projections. 

A section of E' = TQ/SO( 3) = TR 2 x R 3 — > R 2 is a pair (X,u), where A is a 
vector field on R 2 and u : R 2 —> R 3 is a smooth map. Therefore, a global basis of 
sections of TR 2 x R 3 —> R 2 is 

= eS = ( l’ 0) ' 

e 3 = (0, ui), e' 4 = (0, u 2 ), eg = (0 ,u 3 ), 


where u±, u 2 , u 3 : R 2 —* R 3 are the constant maps 


ui{x,y) = (1,0,0), u 2 {x,y) = (0,1,0), u 3 (x, y) = (0,0,1). 


In other words, there exists a one-to-one correspondence between the space Sec (E' = 
TQ/SO( 3)) and the G-invariant vector fields on Q. Under this bijection, the sec- 

d d 

tions e\ and e( correspond with the vector fields — and — and the sections ei, e\ 

ox ay 

d d d 

and e'r correspond with the vertical vector fields -r—, —— and ——, respectively. 

oqi oq 2 dq 3 

The anchor map p’ : E' = TQ /50(3) = TR 2 xl 3 ^ TR 2 is the projection over 
the first factor and, if [-, ■]' is the Lie bracket on the space Sec (E' = TQ /50(3)), 
then the only non-zero fundamental Lie brackets are 


[4 1 


41 / = e 5. 14)41' = e 3> 


[4.41' = 


= e 4 . 


From (17.91) and (17.101) , it follows that the Lagrangian function L = T and the con¬ 
straint submanifold M are 50(3)-invariant. Consequently, L induces a Lagrangian 
function L' on E' = TQ/SO( 3) and, since M is closed on TQ, the set of orbits 
M' = M/SO(3) is a submanifold of E' = TQ/SO(3) in such a way that the 
canonical projection <h| m ■ -M. —> M' = M/SO{3) is a surjective submersion. 

Under the identification between E' = TQ/ 50(3) and TR 2 x R 3 , L' is given by 


L'(x,y,x,y;u}i,W2,U3) 
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where (x,y,x,y) and (x>i, u 2 , W 3 ) are the standard coordinates on TR 2 and R 3 , 
respectively. Moreover, the equations defining M' as a submanifold of TR 2 x R 3 
are 

x — rui 2 + fly = 0, y + rwi — Six = 0. 

So, we have the constrained Lagrangian system (L', M.') on the Atiyah algebroid 
E' = TQ/ 50(3) = TR 2 x R 3 . Note that the constraints are not linear, and that, 
if A' is the Liouville section of the prolongation T E E' , then the restriction to M.' 
of the vector field (p') 1 (A') is not tangent to M’. 

Now, it is clear that the tangent bundle TQ = TR 2 x T(SO{ 3)) = TR 2 x 
(S'0(3) x R 3 ) is the total space of a trivial principal S'0(3)-bundle over E' = 
TQ / SO{3) = TR 2 x R 3 and, in addition (see [13! Theorem 9.1]), the prolongation 
T e 'E' is isomorphic to the Atiyah algebroid associated with this principal SO( 3)- 
bundle. Therefore, the sections of the prolongation T E E' —> E' may be identified 
with the £0(3)-invariant vector fields on TQ = TR 2 x (SO( 3) x R 3 ). Under this 
identification, the constrained dynamics Tfor the system (L/M 1 ) is just 
the 5'0(3)-invariant vector field = P{Tl\m)- We recall that if : TQ —> 

TQ/SO( 3) is the canonical projection, then 

^°r M) = r (i ,^o$. ( 7 . 11 ) 

Next, we give a local description of the vector field (p') 1 (T(u,M')) 011 E' = 
TQ/ 50(3) = TR 2 x R 3 and the nonholonomic bracket {•, -Y nh for the constrained 
system Ad'). For this purpose, we consider a suitable system of local coordinates 
on TQ / SO(3) = TR 2 x R 3 . If we set 


x' = x, y' = y, 

7 r( = rx + k 2 u 2 , tt ' 2 = ry — k 2 uji , ir' 3 = k 2 co 3, 

7 N = ( fc 2+ r 2) (x - ru 2 + fly), = (fc2 + r 2) (y + rwi - fix), 

then (x', j/', 7 r(, 7 T 2 , 7 Tg, 754 , Tg) is a system of local coordinates onTQ/SO(3) = TR 2 x 
R 3 . In these coordinates the equations defining the submanifold AT are n/ = 0 
and 77-5 = 0, and the canonical projection <1> : TQ —> TQ/SO( 3) is given by 

7 x 2 , 773 , 774 , 755 ) = (x,y; 7 Ti, 7 t 2 , 7 r 3 , 7 r 4 , 755 ). (7.12) 


Thus, from (17.71) and (17.111) . it follows that 

(p') 1 ( r u',M')) = (*'^7 + \ m '’ 

or, in the standard coordinates (x, y, x, y; u>i, uj 2 , W 3 ) on TR 2 x R 3 , 


{p ') 1 (f (L 1 ,M’)) — {£( 7 ^ + 

d 

+ y(ir - 


fir 


d 


flk 2 d 
( k 2 + r 2 ) dy ( k 2 + r 2 ) dui± 

flk 2 d fir d 


) 


)}| 


M' • 


y dy ( k 2 + r 2 ) dx ( k 2 + r 2 ) duj 2 

On the other hand, from (17.81) . (17.121) and Theorem 17.91 we deduce that the only 
non-zero fundamental nonholonomic brackets for the system ( L',A4') are 


W^'i) nh= r ’ 

WliAYnh = ^ 

Ynh= ( k 2 k +r 2f 2 

Therefore, we have that 


rk 2 fl , 
(k 2 +r 2 ) X ‘ 


{y'^'iYnh = u 

p 

{^2^ 3 U- (A .2 +r 2)^l 


/' = (^(^(/O + i/'^'lU, for /' € C°°(M'), 


rk 2 fl , 
(, k 2 +r 2 ) V ’ 
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where (p') 1 (Rl') is the vector field on M' given by 


, Al/ _ x r . . d . d , 

^ ^ L '^ ~ ^ 2 +r 2 ^ ^ + (fc 2 




+ z/ , ( 7r 3 - fc2 ^)^7 - + ’’V)^)}! 


(?7T, 


X' 


8. Conclusions and outlook 

We have developed a geometrical description of nonholonomic mechanical sys¬ 
tems in the context of Lie algebroids. This formalism is the natural extension of 
the standard treatment on the tangent bundle of the configuration space. The 
proposed approach also allows to deal with nonholonomic mechanical systems with 
symmetry, and perform the reduction procedure in a unified way. The main results 
obtained in the paper are summarized as follows: 

• we have identified the notion of regularity of a nonholonomic mechanical 
system with linear constraints on a Lie algebroid, and we have characterized 
it in geometrical terms; 

• we have obtained the constrained dynamics by projecting the unconstrained 
one using two different decompositions of the prolongation of the Lie alge¬ 
broid along the constraint subbundle; 

• we have developed a reduction procedure by stages and applied it to non¬ 
holonomic mechanical systems with symmetry. These results have allowed 
us to get new insights in the technique of quasicoordinates; 

• we have defined the operation of nonholonomic bracket to measure the 
evolution of observables along the solutions of the system; 

• we have examined the setup of nonlinearly constrained systems; 

• we have illustrated the main results of the paper in several examples. 

Current and future directions of research include the in-depth study of the reduction 
procedure following the steps of Eli for the standard case; the synthesis of so- 
called nonholonomic integrators mmm for systems evolving on Lie algebroids, 
and the development of a comprehensive treatment of classical field theories within 
the Lie algebroid formalism following the ideas by E. Martinez [M] , 
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